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' A large set of techniques needed to compute decay rates at the two-loop level are derived and 

systematized. The main emphasis of the paper is on the two Standard Model decays 77 — > 77 
and H — > gg. The techniques, however, have a much wider range of application: they give 
practical examples of general rules for two-loop renormalization; they introduce simple recipes 
Q!^ I for handling internal unstable particles in two-loop processes; they illustrate simple procedures 

for the extraction of coUinear logarithms from the amplitude. The latter is particularly relevant to 
I show cancellations, e.g. cancellation of coUinear divergencies. Furthermore, the paper deals with 

the proper treatment of non-enhanced two-loop QCD and electroweak contributions to different 
physical (pseudo-)observables, showing how they can be transformed in a way that allows for 

■ a stable numerical integration. Numerical results for the two-loop percentage corrections to 
H I if — > 77, gg are presented and discussed. When applied to the process pp-^gg + X^H + X, 

■ the results show that the electroweak sealing factor for the cross section is between —4% and 
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1 Introduction 

In this paper we have coUected and systematized a large set of techniques needed to compute decay rates 
and production cross sections at the two-loop level in a spontaneously broken theory characterized by a large 
number of scales. Although the main emphasis will be on the Standard Model (SM) Higgs-boson decays to 
two photons and two gluons and on the related production processes, our techniques have a broader range of 
application. Firstly, they represent practical applications of general rules for two-loop renormalization and 
introduce simple recipes for handling internal unstable particles. Secondly, they illustrate useful procedures 
for the analytical extraction of collinear logarithms: we explicitly prove that the complete H amplitude 
is free from collinear logarithms in spite of the fact that single Feynman diagrams show a collinear-divergent 
behavior. Finally, our techniques allow to represent the non-enhanced two-loop QCD and electroweak con- 
tributions to different physical (pseudo-)observables in a form suitable for a stable numerical integration. 

The lastly mentioned feature represents the important and difficult part of our computation. Since we 
do not rely on any kind of approximation, we can safely cover all kinematical regions including different 
normal thresholds, like the WW one, which are peculiar of Higgs- and vector-boson decays. Therefore, we 
are able to place the electroweak component on the same footing as the QCD one: a non-trivial but essential 
task, since in several situations the size of the pure electroweak higher order corrections is comparable to 
the theoretical uncertainty related to the parton distribution functions. 

Note that two-loop electroweak corrections to 1 ^ 2 processes are currently investigated by several 
groups: the authors of Ref. ^1 have derived the dominant contributions to the H bb decay width, and 
two different collaborations have computed the fermionic [2, and bosonic |3l4j corrections to the effective 
electroweak mixing angle. 

With the advent of the LHC era, it is clear that the production process gg ^ H and the decay mode 
H — !■ 77 are going to play a crucial role for a precise comparison of the experimental data with the SM 
predictions. 

Gluon fusion is the main production channel of the SM Higgs boson at hadron colliders, and both 
virtual and real corrections have been thoroughly investigated since the beginning of the 90's. The QCD 
next-to-leading order (NLO) radiative corrections to the total Higgs-production cross section have been 
first computed below the ti threshold in Ref. [5 , and using an effective-Lagrangian approach, where the 
top quark is integrated out, in Ref. 6J. Formally, the second method defines the heavy-top limit, and it 
has represented the starting point for subsequent higher order improvements. A one-dimensional integral 
representation which proves well-set for numerical evaluation has been later derived in Refs. |7l8j for the 
entire Higgs-mass range; moreover, the virtual component has been evaluated analytically in a closed form 
by the authors of Ref. [9] . 

Recently, two different groups |10|llj have provided independent checks of the work of Ref. f9] , general- 
izing the analytic result in a framework where the coupling of the Higgs to the external particles is mediated 
by a scalar field (see also Ref. [12] )■ Since QCD NLO corrections increase the cross section by more than 
70 %, there was a flurry of activity on higher order QCD effects. Partial results for the next-to-next-to- 
leading order (NNLO) corrections to the total cross section in the heavy-top limit have been obtained in 
Ref. [13] for the two-loop virtual corrections to the effective heavy-top Hgg vertex, and in Ref. [14] for the 
soft components; the complete NNLO result has been finally derived by three different groups in Ref. [TS] . 
and later supplemented by an all-order resummed calculation of multiple soft-gluon emission at next-to- 
next-to-leading logarithmic (NNLL) accuracy in Ref. [H] and including the fourth logarithmic (N^^LL) order 
in Ref. TT] (see also Ref. [H]). In addition, the effect of a jet veto on the inclusive cross section has been 
studied at NNLO in Ref. [19j . and an improvement with respect to the heavy-top limit has been recently 
obtained by the authors of Ref. [2^. Because of the cuts on the final states typical for hadron-coUider 
phenomenology, fully differential perturbative results have a privileged role; the differential cross section for 
Higgs production has been derived through NNLO in QCD in Ref. [211, and later cross-checked through an 
independent subtraction formalism by the authors of Ref. [22] . Note that the techniques used in Ref. [2T] 
have been recently extended in Ref. [23] to compute the QCD NNLO cross section for the H WW Ivlv 
signal at the LHC. The Higgs-production mechanism in the gluon- fusion channel has clearly a close affinity 
with the Higgs decay to two gluons; the state of the art in QCD is presently represented by the two works 
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of Ref. [21], with the fuU NNLO resuh and the N^^LO calculation in the heavy-top limit respectively. 

While QCD corrections to Higgs production through gluon fusion at hadron colliders are well under 
control, electroweak effects are less understood. The NLO corrections to the total cross section were evaluated 
in Ref. [25l26j in the heavy-top limit, and turned out to be less than 1%. The contribution due to the 
light fermions has been calculated analytically in Ref. (27], and found to be more sizable; the remaining 
component of the amplitude involving the top quark has been computed by means of a Taylor expansion 
in the Higgs external momentum in Ref. [28]. It is worth stressing that the work of Refs. [27|28j has not 
been independently cross-checked by a different group; moreover, for obvious reasons, the validity of the 
result of Ref. [28j is restricted to the kinematical region below the WW threshold. Under the hypothesis of 
factorization with respect to the dominant QCD soft and coUinear radiation, the impact of the electroweak 
corrections to Higgs-production in proton-proton collisions has been estimated in Ref. |29j . 

The decay of the SM Higgs boson into two photons has also a great phenomenological interest: on the 
one hand, it provides a precious information for the discovery at hadron colliders [3^; on the other hand, 
an upgrade option at the planned ILC will allow for a high-precision measurement of the partial width into 
two photons [31], leading to quantitative tests on the existence of new charged particles. QCD corrections 
to the decay width of an intermediate-mass Higgs boson have been computed at NLO in Ref. [32] and at 
NNLO in Ref. [33]; the NLO result has later been extended in Ref. [34] to the whole Higgs-mass range. 

Electroweak NLO corrections have been computed in Ref. [35|26j in the large top-mass scenario (the 
so-called "dominant" corrections) and in Ref. [3S] in the large Higgs-mass scenario. The complexity of this 
complicated (multi-scale) computation is reflected into the fact that seldom an a priori dominance lives 
up to its promise; for instance, O {GpMf) corrections in the range 80 — 150 GeV do not dominate at all. 
Recently, the two-loop contributions induced by light fermions have been derived in Ref. [57] , and electroweak 
corrections due to gauge bosons and the top quark have been given in term of an expansion in the Higgs 
external momentum in Ref. [37] , a result subsequently extended in Ref. [38] to cover all kinematical regions 
including the notoriously difficult WW threshold and the NLO QCD corrections. 

In this paper we show details about the evaluation carried on for the H —> 77 decay in Ref. j38j and 
generalize the result to the gg ^ H process to the entire Higgs-mass range. Since we are not bound to rely 
on any kind of expansion, neither in the bosonic sector nor in the top-bottom one, we can produce extremely 
accurate results for any value of the Higgs-boson mass, including the full dependence on the W-, Z- and 
Higgs-boson masses and on the top-quark mass. A consistent and gauge-invariant treatment of unstable 
particles allows to produce precise results around the WW threshold. 

In our approach, Feynman diagrams, up to two loops and including QCD, are generated by means of an 
automated procedure which does not rely on any external package, and is implemented in the FORM |3!?j 
program QraphShot [ID]. After projecting the amplitude for a given (pseudo-) observable onto form factors 
and taking traces over Dirac matrices, three basic simplifications are done recursively: at first, reducible scalar 
products are removed; next, tadpole integrals are reduced using integration-by-parts identities (IBPIs) [41^; 
finally, the symmetries of each diagram are exploited in order to reduce the number of integrals to be 
evaluated. 

The second step in the calculation of any two-loop process concerns renormalization; here we heavily rely 
on the results of Refs. |42|43|44] . After assembling diagrams we perform the usual canonical tests to ensure 
that our result is correct, checking all possible unrenormalized Ward-Slavnov- Taylor (WST) identities. The 
next logical steps consist in removing all ultraviolet divergencies, with the due caution to the problem of 
overlapping divergencies, and checking the renormalized WST identities. Then, we express renormalized 
parameters in terms of physical observables belonging to some input-parameter set. 

The third step consists in the analytical extraction of coUinear logarithms of Feynman diagrams. It is 
worth noting that the amplitude for ^ 77 is coUinear free, and one could adopt the approach where all 
light fermions are considered massless; therefore, the coUinear behavior of single components is controlled 
in dimensional regularization and coUinear poles cancel in the total. We prefer another approach, where 
coUinear singularities are controlled by the light-fermion masses; partial components of the complete result 
are divergent, but we check that all logarithms of coUinear origin cancel in the complete answer. 

Finally, the remaining coUinear-free contributions are written in terms of smooth integral representations 
using the methods developed in Refs. |45I46I47I48I49| . and then evaluated numerically. 
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The outline of the paper is as follows: in section [5] we summarize our notation and conventions. Details 
about the construction of the H —f "/"/{gg) amplitude, including renormalization and choice of the input- 
parameter set, are discussed in section [3] In section 2] we explain how we manipulate Feynman integrals 
in order to simplify the amplitude, removing scalar products, reducing tadpole integrals and symmetrizing 
Feynman diagrams. The behavior of two-loop diagrams around a normal threshold is thoroughly investigated 
in section [5l In section [6] we discuss the extraction of coUinear singularities, and in section [7] we describe 
details about the evaluation of massive diagrams. Numerical results are shown and discussed in section [8l 
Finally, we summarize our conclusions in section [9l Additional material is contained in the appendices. 

2 Notation and conventions 

Regularization. We employ dimensional regularization [50| . denoting the number of space-time dimensions 
hy n — A — e, and we define short-hand notations for regularization-dependent factors, 

A„v = 7 + liiTT + In ^5", ^uv{x) = A[jv - In , (1) 

where 7 = 0.5772156 • • • is the Euler constant, ^ is the 't Hooft mass unit, M stands for the bare or 
renormalized VF-boson mass (in the following we will not distinguish unless strictly needed) and a; is a 
positive-definite kinematical variable. In our conventions the logarithm has a cut along the negative real 
axis and it is understood that for all masses ^ — iQ. 

The structure of poles in dimensional regularization is better exploited in terms of universal factors. 
Any one-loop integral can be formally decomposed as ( {^} ) = X]fe=-i ( {'} ! ^ ) F^ix), where {1} 
stands for a set of arguments representing powers of inverse propagators, external kinematical variables and 
masses of internal particles and x denotes a scale which depends on the kind of integral: M for tadpoles 
and a squared external momentum for other configurations. The dependence on the dimensional regulator 
e and the regularization-dependent factors introduced in Eq.([T]) is encapsulated in the universal one- loop 
ultraviolet (hereafter UV) factors 

F\{x) = -^-^AMx) + ^Al,{x)e, F^{x) = 1 - ^ A^^x) e, F^{x)^e. (2) 

It is worth noting that, because of overlapping divergencies due to UV-divergent one-loop sub-diagrams, we 
include O (e) terms in all one-loop results. 

A generic two-loop integral can be written as ({0) — J2k=-2 ( {0 ! ^) F^{x), with two-loop 
UV factors given by 

Fl,{x) = ^-^^ + ^AlAx), Fl,{x)^^-A„Axh F,'{x)^l. (3) 

Note that the product of two one-loop integrals can be written in terms of the same UV decomposition of a 
genuine two-loop integral. 

Finally, let us define our soft/coUinear factors; for the electroweak part of the calculation they are equal 
to the UV factors of Eq.([3]). In appendix [Cl we will extend their definition to cover non-abelian QCD 
configurations. 

Classification of two-loop integrals. In this paper we frequently use the notion of scalar, vector and tensor 
two-loop Feynman integrals. An arbitrary two-loop scalar diagram can be cast as 

G'^'-'-^ d-q,d-q,Y[ik^+m^)-' Y[ (k^ + mj)-' 11 i^f + mfr\ (4) 

i—l j—a+1 l—a+j+1 

where a,/3 and 7 give the number of internal lines containing the integration momenta (71,(72 and qi — (72, 
and we have introduced 

ki^qi + Vh Pj' i = 1, . . . , a, 

h = qi-q2+J2j=iV}^P3^ i = a + l,...,a + 'y, (5) 
h^q2 + EjLi vfjPj, i^a + j + l,...,a + j + f]. 
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Here N is the number of vertices connected to external lines, r/fj = ±1,0 and {p} is the set of external 
momenta. The capital letter G identifies the number of external legs: T for tadpoles, S for self-energies and 
V for vertices. 

The triplet of numbers aPj can be represented in an extremely compact way: first, we introduce k = 
7max [amax (/? " 1) + « " 1] + 7; where for tadpoles amax = 1 and 7inax = 1, for self-energies amax = 2 and 
7max = 1 and for vertices amax = 2 and 7max = 2; next, we associate a letter of the alphabet to each value 
of K and we derive a set of alpha-numerical correspondences, 



G = T lll^yl, G = S => 111^ A, 121 ^ C, 131 ^E, 221- 

G = V => 121^ E, 131 ^ /, 141 ^ M, 221 ^ G, 231 ^ K, 222 H. 



(6) 



All the different groups (sometimes called families or topologies) of Feynman diagrams with different num- 
bers of external and internal legs have been classified for completeness in Figs. [551 HH EZl and [55] of ap- 
pendix]^ where the graphical relation between the diagrams and the corresponding integrals can be found. 
Note that we follow the notation of Ref. 51j for one-loop integrals. At the two-loop level we identify 
one tadpole (vacuum) topology, T'*, four self-energy topologies, S^, S"^, S'^, S° , and six vertex topologies, 
V^, V , y", V^, , V" . Note that factorized two-loop topologies, associated with the product of two 
one-loop Feynman integrals, do not receive a particular name. 

For a detailed analysis of scalar two-loop self-energies and vertices we refer the interested reader to 
Refs. [45l46l47l48j . The presence of non-trivial structures containing integration momenta in the numerators 
of two-loop integrals requires to introduce tensor structures and form factors, as described in sections 7 and 9 
of Ref. [49] for higher rank self-energies and vertices. 

Miscellanea. In order to keep our results as compact as possible, we introduce a short-hand notation for 
integrals over a simplex of Feynman parameters. 



/ 



dSn{{x})f{xi,- ■ •,a;„) = Jl 
where xq 



dxif{xi,---,x„), / dC„{{x})f{xi,---,x„)^ J|da;i/(a;i, 



1 " 



(7) 



1. In addition, the so-called ' 



dx- 



dx 



x ^ a 
f{x,{z}) In" a; 



and ' 



= / dx 



\-' distributions will be extensively used, 
fix,{z})^fia,{z}) 



= dx 



X — a 

[f{x,{z})-f{0,{z})] In-x 



a = 0,1, 



''dxdy^^^^y^^^^^ 
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xy 



= dxdy 



fix, y, {z}) - fix, 0, {z}) - /(O, y, {z}) + /(0, 0, {z}) 
xy 



(8) 
(9) 



3 The amplitude for H ll{gg) 

We consider the decay of an arbitrary Higgs boson into two photons h (— P) -I- 7 (pi) -I- 7 (^2) 0, where 
P = Pi + P2 (all momenta are incoming) . The amplitude A can be written as 

A = Z^^ el ^^,(0, 0, -Ml), (10) 

where = e^ipi,Xi), with i = 1,2 and = ±, are the photon polarization vectors and .E^ and Z„ are the 
photon and Higgs-boson wave-function renormalization factors. Afj_i, is the amputated Green's function for 
h 77, whose tensor structure reads 



9 

G - 
16772 



A,ApIpIP^) = ^ FM,pIP^)5^.+ J2 F^''\plplP^)p.,^p,, + F,iplplP')e^.^0p'}p^. (11) 
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li h ^ H, the Standard Model CP-even Higgs boson, then — and G — Sg, where g is the bare 
S'[/(2) -couphng constant and sg — sin 6* {cg = cos 6*) is the sine (cosine) of the bare weak-mixing angle. The 
form factors Fj^, Fp^''' and are functions of the off-shell kinematical invariants pi, p\ and P^; they have 
to be extracted with suitable projection operators and can be subsequently evaluated for on-shell external 
momenta, p\ = P2 — Q and — —M"^, where is the on-shell Higgs-boson mass. A general framework 
for projection operators is discussed in appendix [Bl following the approach of Ref. [52j . 

Note that the amplitude for the decay of the SM Higgs boson into two gluons, H gg, is obtained from 
Eq. (|10p and Eq. (|lip . replacing Za with Zg, the gluon wave- function renormalization factor, and introducing 
the appropriate color indices for the gluon fields; in addition, the overall factor g^ Sg of Eg. pTjl is replaced 
by g g^, where gs is the S'f7(3) -coupling constant of strong interactions. 

Because of the absence of a tree-level H^^{gg) coupling in the SM, the lowest-order amplitude is generated 
by one-loop fermionic and bosonic diagrams, as shown in Fig.[T] note that the virtual particles do not decouple 
also for masses larger than the Higgs one, since their coupling to the Higgs boson grows with their mass. 




Figure 1: Examples of one-loop diagrams for the decay H — > 77. Other bosonic diagrams, not shown here, contain 
charged unphysical Higgs-Kibble scalar and Faddeev-Popov ghost fields. For H gg, only fermionic diagrams 
appear. 

Having singled out the sg dependence due to the couplings of the external photons, the one-loop form 
factors do not show any residual dependence on 9, because one-loop diagrams containing Z bosons or neutral 
unphysical Higgs-Kibble scalar and Faddeev-Popov ghost fields do not show up. Therefore, form factors at 
one loop contain only the bare masses of the W boson, m„, the Higgs boson, m^, and the fermions, nif. 
This simple observation will allow in the following for a straightforward implementation of two different 
renormalization schemes. On the one hand, the on-shell masses of the W boson, M^^., the Higgs boson, M^, 
and the fermions, Alf, will be part of any renormalization scheme, and the associated derivatives of the one- 
loop form factors will have to be evaluated. On the other hand, g and sg, collected as simple pre-factors in 
Eq. (jll[) . will be connected to Gp, the Fermi-coupling constant, and a, the fine-structure constant, or M^, the 
on-shell Z-boson mass. The choice of the two different input-parameter sets {Gp, a) or {Gp, M^) will define 
our renormalization scheme. Finally, the 5L/(3)-coupling constant gs will be related to the strong-coupling 
constant as(/i^) evaluated at the appropriate renormalization scale /i^. 

The form factor is present at two loops in diagrams containing a fermion sub-loop, as shown in Fig. [21 
but it does not arise at one loop due to the lacking of axial fermion couplings. Due to CP invariance the 
contribution of F^ vanishes in the total. Therefore, we can circumvent the notorious problems associated 
with the definition of the 75 matrix in dimensional regularization [50j (see a detailed discussion in Ref. |53j ) 
employing a completely anticommuting 75. 




w ^ t "f w 



Figure 2: Representative two- loop diagrams giving a non- vanishing contribution to the form factor F^, because of the 
presence of an axial coupling between a massive vector boson and a virtual fermion pair. 

Finally, we observe that Bose symmetry sets four constraints on the other form factors for off-shell 
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external particles, 

{pIpIp') = {pIpIp') , {pIpIp') = {pIpIp') , 

{pIpIp') = 4^^) {pIpIp') , Fj,''^ {pIpIp') = F^'^^ (P^.P?,/'^) . (12) 

Although Fp , Fp and Fp ' are irrelevant for the calculation of the amplitude due to the condition 
■ Pi — for on-shell photons (gluons), they play an essential role in computing the relevant Ward-Slavnov- 
Taylor identities (hereafter WSTIs) for H — > 77(53)- 

Before discussing WSTIs in more detail, let us define what we use for the Higgs-boson mass and wave- 
function renormalization factor. 



3.1 Higgs-boson mass and wave-function renormalization 

We evaluate the amplitude for P^ — — A/^, where is the on-shell Higgs-boson mass, defined as the 
location of the pole of the real part of the associated Dyson-resummcd propagator A^, 

Re [A^ip'')]'' \p2^_Ml =0 => Ml = ml - ReS„(-M|). (13) 

Here is the bare mass and Y,h is the sum of all one-particle irreducible (IPI) diagrams for the Higgs 
self-energy. For definiteness, Af^ would be the experimental mass extracted from a fit of the hypothetical 
Higgs-resonance line shape obtained by means of a Breit-Wigner function with a running width, 

^"^^'^ " p^ + Ml+\p^T^/M^- (1^) 

It is well-known that the on-shell definition of mass for an unstable particle is gauge-parameter dependent 
beyond one loop, and it should be replaced by the complex pole, s„ in our case, including real and imaginary 
parts of , 

[A„(/)]"^ |p2=_,^ ==0 ^ s„ = ml - S«(-s„). (15) 

As it is often stated |54j . the complex pole is a property of the S matrix, and therefore gauge-parameter 
independent at all orders in perturbation theory; the proof of this property relies on the Nielsen identities and 
can be found in Ref. [55j . Therefore, fixed-order computations beyond one loop or resummation- improved 
one-loop predictions require to perform mass renormalization dropping the concept of on-shell masses and 
employing complex poles; for a thorough discussion at one and two loops we refer to Ref. |44| . 

Concerning H — > 77 at two loops, the absence of a tree-level amplitude implies that Higgs-mass renor- 
malization has to be performed at one loop: this suggests that an on-shell definition of the Higgs-boson 
mass, being gauge-parameter independent at one loop, proves adequate for our purposes. Nevertheless, we 
will show in the following that a naive implementation of on-shell mass renormalization for the Higgs boson 
breaks the standard two-loop WSTI for 7J ^ 77 above the VFVF-production threshold (light-fermion masses 
are neglected). On the one hand, the bare identity develops an imaginary part above the WW threshold; 
on the other hand, the Higgs self-energy appearing in Eq. (fT5)) . E^, develops an imaginary part at one loop 
for > 2M^, but on-shell mass renormalization selects its real part only. This introduces a mismatch 
between real and imaginary parts at the level of the renormalized WSTI. Therefore, the usual statement that 
renormalization schemes with complex poles or on-shell masses are equivalent at one loop should be taken 
with some caution; for H 77 at two loops, the two schemes are equivalent only below the IFW^-production 
threshold. 

A second delicate point concerns the use of a Higgs-boson wave-function renormalization (hereafter 
WFR) factor introduced a la Lehmann-Symanzik-Zimmermann (LSZ) [56] . The LSZ formalism, in fact, is 
unambiguously defined for stable particles, but it requires some care when external unstable particles appear. 

Let us consider first the case of the photon in QED. Here, the LSZ WFR factor is fixed by the condition 
that the transverse part of the photon Dyson-resummed propagator, A^, dressed by its WFR factor Za, has 
unity residue at the pole in the p'^ plane, 



A2p') |p2=o ^ 2, = 1 - l],,p(0), I].,p(p^) = (16) 



_d_ 
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where is the sum of all IPI diagrams for the transverse part of the photon self-energy (we will show 
in the following that the presence of a photon/^-boson mixing in the SM does not change the structure of 
Eq. (fTB| V Note that it has been shown long ago [57l58j that the LSZ definition of WFR factors leads to 
gauge-parameter independent 5-matrix elements and respects unitarity. 

For a stable particle, the pole of the Dyson-resummed propagator is real; being the Higgs boson an 
unstable particle, the pole of its Dyson-resummed propagator in the plane is complex, and the Higgs 
WFR factor is not unambiguously defined. We define it requiring that the real part of the Higgs-boson 
Dyson-resummed propagator, dressed by its WFR factor Z^, has unity residue at the on-shell Higgs-boson 
mass, 

Z«=ReA [A„(p2)]-i|p,_^,^ ^ Z^ = l- RcS«,p(-M2), ^ (17) 

Other definitions, involving the real part of the derivative of the propagator at the complex pole [59], or the 
complete derivative, including imaginary parts [60j . do not look appropriate for our computation, which deals 
with the ideal situation of an asymptotic state containing the Higgs boson, whose momentum is required to 
be real. 

Note, however, that the definition of a real WFR factor at the on-shell mass is a possible source of 
inconsistencies. The authors of Ref. [61| . for example, have shown that unitarity-breaking terms can show 
up in the final form for the amplitude. From a formal perspective, indeed, transition amplitudes have to 
be defined in terms of asymptotic states containing stable particles only. To this respect, the necessary 
modifications to the theory have been described in Ref. [62j . where unitarity and causality of the resulting 
S matrix are proven. 

In the following, we will use a compromise. At first, we will show that, concerning H — + "f"f{gg) at 
two loops, a naive use of Eq. (fT7|) leads to unphysical singularities at thresholds. Then, we will modify our 
renormalization scheme introducing complex poles. 



3.2 Ward- Slavnov- Taylor identities 

Ward-Slavnov- Taylor identities (WSTIs) [63] for H — > jjigg) are essential for organizing our compu- 
tation. They allow to perform severe checks on the algebraic structure of the amplitude, and set strong 
constraints on the number of independent form factors. Although in the following we will focus on the 
iJ ^ 77 process, it is evident that analogous considerations can be applied to H ^ gg. 

Let us consider the simplest case, the so-called doubly contracted WSTI with an on-shell Higgs boson 
and both off-shell photons. After replacing in Eg.fTO)) the photon polarization vectors with their associated 
four- moment a, we get 



-1/2 u 



pt pi A^M.pi.-K 



32 

1 
2 



'-±Z-"'[[Ml +pI+ pI) [f^^pI-MD + Y: pIF^^'^ {pIpI -Ml) 

i=i 

{Ml +pi+pl) F^^'\pIpI-MI)\ ~ 2pIpIF^''\pIpI-MI)} = 0. (18) 



-1/2 



Here the Higgs boson is assumed on its mass shell, — —Ml, and provided with its WFR factor Zh' 
whereas the form factors are evaluated for off-shell photons and photon WFR factors are consistently not 
included. After showing that all form factors are regular for on-shell photons, one can set Pi = P2 ^ ^^'^ 
get the well-known constraint 



WXde(0,0) 



g^4 ^-1/2 2 



F„(0,0,-Af2) 



^4^^)(0,0,-M^) 



0, 



(19) 



which shows that the H —^ 77 amplitude can be described through a single form factor. Before attempting 

(21) 

any evaluation of the amplitude, we extract both F^, and Fp employing projection operators for off-shell 
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photon momenta. Next, we check that form factors are regular for on-shell photons, and we compute the 
WSTI, verifying that Eq. fTO)) holds at the algebraic level. This provides a strong test on several procedures 
used throughout the calculation of the amplitude, like the reduction of tensor integrals to scalar ones. Note, 
however, that the WSTI, being satisfied at the algebraic level, does not provide any information on the 
analytic structure of the loop integrals themselves. Here, a severe test on the result will be supplemented by 
the extraction of all coUinear logarithms: they have to cancel in the full amplitude, and indeed we will show 
that they cancel. 

Two additional relations among the form factors can be readily obtained by considering simply contracted 
WSTIs, where one off-shell photon leg is contracted with its four-momentum, and the other one is saturated 
with the polarization vector, put on the mass shell and provided with its WFR factor. 
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16 7r2 
1 



/ N -7-1/2 -7-1/2 

(pi • 62 ) Zh 



iMl+pi)Fi''\pi,0^~Ml) 



FApi,0, -Ml) + pi Fi''\pl 0, -Ml 



(20) 



WIscMp': 



-1/2 ^-1/2 fj. 



9^ 4 
16 7r2 



/ N -7-1/2 -7-1/2 

(P2 • ei) Za ' Zh ' 



FUO,pI-MI) + piFi''\0,pl-Ml) 



-i^{mI+pI)f, 



(21) 



= 0. 



(21) 



In Fig. [3] we show diagrammatically Eq. p8)) , Eq. (pO)) and Eq. ((2T|) . In addition, we introduce the notions 
of: 1) physical source: the external on-shell leg is multiplied by the appropriate WFR factor, and, for the 
case of photons, supplemented by the polarization vector; 2) contracted source: the external off-shell photon 
leg is contracted with its four-momentum. 



H 



y ~yww# 



= 



H 



Figure 3: Diagrammatic representation of Ea. (|18p (doubly contracted WST identity), Ea. (|20p and Ea. (|2ip (simply 
contracted WST identities). Physical sources are denoted by a gray box, contracted ones, for photons, by a black 
circle. In the latter case the on-shell relation pf = is relaxed. 



The computation of a doubly or simply contracted WSTI for an on-shell Higgs boson beyond the leading 
order approximation requires to pay special attention to the interplay between the bare and on-shell Higgs- 
boson masses. Let us consider the doubly contracted WSTI of Eq. (fT9| at two loops. The full set of diagrams 
can be organized in three classes: 1) two-loop diagrams where the Higgs-boson WFR factor is set to unity, 
Zh — 1, and the bare Higgs-boson mass is identified with its on-shell experimental value, m„ — Mh\ 2) 
one- loop diagrams supplemented with Zh evaluated at one loop by means of Ea. (|17p . where we employ 
again the tree-level relation mn — Mh] 3) one-loop diagrams where Zh = 1, and we replace the tree- level 
identity between the bare and the on-shell Higgs-boson masses with the solution of the Higgs-boson mass 
renormalization equation of Eq. ()13p . with the Higgs-boson self-energy Y]h evaluated at one loop. 

In Fig. m we show a representative diagram for each class; let us consider in detail their expansions in 
terms of the bare 5'C/(2)-coupling constant g. The first diagram is 0{g^); the second one is 0{g^), and gives 
an 0{g^) contribution once provided with the WFR factor Zh at one loop. The most subtle point is related 
to the third diagram, which is 0{g^) and is multiplied by a tree-level WFR factor, Zh = 1. However, it 
entails a vertex where one Higgs boson and two unphysical charged Higgs-Kibble scalar fields are coupled 
proportionally to m^. Since we put the Higgs boson on its mass shell, the bare mass has to be removed, 
introducing everywhere the on-shell value Mh- Therefore, first we use Eq. p^ and evaluate T,h at order 
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0{g^); next, we replace the solution in the third diagram and obtain an 0{g^) contribution, which proves 
essential in fulfiUing the WSTI. 



W 7 




WFR T 



7 

D-----<^ 



Figure 4: Representative diagrams for the doubly contracted WSTI of Ea. (|18|) at next-to- leading order. As in Fig. [3] 
black circles correspond to off-shell photons contracted with their four-momenta, and a gray box stands for a physical 
Higgs-boson source. Inclusion of the Higgs-boson WFR factor Zh in the second diagram and Higgs-boson mass 
renormalization in the third one are both performed at one loop. 

Finally, we stress that both doubly and simply contracted WSTIs assume that the Higgs boson is emitted 
by a physical source. However, as shown by 't Hooft and Veltman in their seminal paper |58j . a broader 
class of WSTIs can be obtained when we inspect, at the diagrammatic level, the effect of an infinitesimal 
shift of the gauge-fixing functions, C — > C -I- e 7?,, with e ^ 0. Here a denotes one of the gauge fields 
and 7?, stands for an off-shell source emitting one or more off-shell particles. As shown in Ref. [58 , off- 
shell WSTIs connect physical or contracted sources to special sources obtained by: 1) subjecting Tl, to an 
infinitesimal local gauge transformation; 2) replacing the parameters of the transformation with the related 
Faddeev-Popov ghost fields. Let us consider a simple example one can derive for H — > 77. After choosing 
the gauge-fixing function for the photon field in 't Hooft-Feynman gauge, = — we introduce 
two off-shell sources emitting one photon or one Higgs boson, TZa = JaA^ and TZh — JhH . Next, we 
perform a local SU{2) x U{1) gauge transformation, and replace the infinitesimal gauge parameters of the 
transformation by the appropriate Faddeev-Popov ghost fields. We derive two special sources, 



TZa ^ J'^ [A^ + igsg {X-W+ - X+W;) - a^F^J, 

n„ ^ j„\h + ^Y"" + ^ {x-ip+ + x+ip-) 



(22) 



where Y"^, Y^, X^ and X^ are the Faddeev-Popov ghost fields and <y9°, (p~^ and ip~ are the unphysical 
Higgs-Kibble scalar fields. Off-shell WSTIs for ^ 77 are then readily obtained connecting at least one 
special source with physical or contracted sources. In Fig. [5] we show two examples: in the first case, the 
Higgs boson and one photon are emitted by physical sources, and the second off-shell photon is connected 
to a special source; in the second case, both off-shell photons are contracted with their four-momenta, and 
the off-shell Higgs boson is emitted by a special source. In this case, WFR factor and mass renormalization 
for the Higgs boson have not to be included. 



3.3 One-loop counterterms 

In this section we shortly summarize three aspects of the renormalization procedure which are needed 
before performing finite renormalization: tadpole renormalization, diagonalization of the neutral sector and 
ultraviolet counterterms. A detailed analysis can be found in Refs. [42|43j . 

We fix the gauge for the electroweak sector of the SM Lagrangian introducing six gauge parameters 
{t = A, Z,AZ, 



-af 



C^C 



— ^ 'T—df^Ap^ — ^Az dp.Zp,, 



1 

2 

ce ?H' 



Here Af^ , and are the fields for the photon and the Z and the W bosons and ip^ and ip^ are the fields 
for the neutral and charged unphysical Higgs-Kibble scalars. For the QCD sector, we employ the following 
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H 




7^7 7 



= 



Figure 5: Off-shell WSTls involving special sources, denoted by gray ovals. The special source for the photon field can 
emit a Faddeev-Popov (FP) ghost field (dot-line) and two W-boson/charged FP field (wavy/dot-lines) couples. The 
special source for the Higgs boson can emit three neutral and charged Higgs-Kibble/FP field couples. See Eq. (|22p for 
the explicit expressions of the special sources. As usual, contracted sources are denoted by black circles and physical 
ones by gray boxes. 

choice, 

E^''^"^'''"' C^''^ = -7^5^G^, (24) 
which for = 1 reduces to the usual 't Hooft-Feynman gauge. 
3. 3. J Higgs tadpoles 

Tadpoles in spontaneously broken theories have been discussed by many authors (see Refs. |64l65j ). The 
tadpole-renormalization prescription used in this paper will be the (3t scheme, thoroughly discussed in sec- 
tion 2.3 of Ref. [Hj. 

Following notation and conventions of Ref. , the minimal Higgs sector of the SM is provided by the 
Lagrangian 

Cs = -iD^K)\D^K) - ^i^K^K - iX/2)iK^Kf, (25) 



where the covariant derivative is given by 



where g' /g — —se/cg, r" are the standard Pauli matrices, is a triplet of vector gauge bosons and a 
singlet. For the theory to be stable we must require A > 0. In addition, we choose /^^ < in order to have 
spontaneous symmetry breaking. The scalar field belongs to the minimal realization of the SM, where C and 
the Higgs-Kibble fields ipo, ^pi and Lp2 are real. In particular, we choose C -I- i(fQ to be the component of K 
to develop the non-zero vacuum expectation value, and we set ((^o)o = and (C)o 7^ 0. We then introduce 
the (physical) Higgs field as H = C, — v. The parameter v is not a new parameter of the model and its value 
must be fixed by the requirement that (i?)o = (i.e. (A')o ~ [1/ \/2){v,Q))^ so that the vacuum does not 
absorb or create Higgs particles. 

It is then convenient to define the bare parameters (the W boson mass), mu (the mass of the physical 
Higgs particle) and (3t (the tadpole constant) according to the following "/3t scheme" , 

V ^2m,y{l+ [3t)/g 

A = (gTO„/2TO„,)^ . (27) 
= + I (2m„V.g)" I ^? = -\ml 
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Next, we define /3t = Pt^ + /3ti3^ + /3t2ff^ + • ■ • and we fix the parameter Pt such, that the vacuum- 
expectation value of the Higgs field remains zero order by order in perturbation theory. At the lowest order 
we can simply set Pt = 0, i.e. j3ta = 0. At one loop we take into account all tadpole diagrams and we get 



= 16^r A. [7^4 + 7< «o(™"-)J + 1. < 

where mz = mw/ce is the bare mass of the Z boson. The tadpole integral ao reads 

1 

ao{m) = ^ ao{m] k)Fl{inn?„), (29) 



with universal one- loop UV factors given by Eq.([2]) and coefficients 



, X 1 r 1 / 1 , w? 



ao(m;-l) = -2, ao(m; 0) = -1 + ln — , aqK 1) = - -1--C(2)+ l--ln— In— , (30) 

2 L 2 V 2 m„,- / m„,- J 



where C,{x) is the Riemann's zeta function. The full list of /3t-dependent Feynman rules needed for our 
computation is given in appendix B of Ref . [12] . 

3.3.2 Diagonalization of the neutral sector 

The Z-j transition in the SM does not vanish at zero squared momentum transfer. Although this fact does 
not pose any serious problem, not even for the renormalization of the electric charge, it is preferable to use 
an alternative strategy. We will follow the treatment of Ref. [Blj; consider the new 5't/(2)-coupling constant 
g, the new mixing angle 9 and the new W mass to„' in the /3t scheme, 

-n ,-r\ ' sine _ 2rn^y fgm„Y ^ 1 2 /oi^ 

cosy g \2m„J 2 

where T = Tig^ + T2g'^ + ■ ■ ■ is a new parameter yet to be specified. This change of parameters entails new 
and fields related to and i?° by 



Zu\ ( cos Q — sin 0\ f B 



A^J cos 9 ) \Bl 



(32) 



In our approach T is fixed, order by order, requiring that the Z-7 transition is zero at = in 't Hooft- 
Feynman gauge. The explicit result for Fi is 



ri = " \ ao{m„), (33) 

16 TT^ 



with ao(m) given in Eq. (|30p . The full list of F-dependent Feynman rules needed for our computation is given 
in appendix C of Ref. [42 . In addition, for iJ ^ 77 we need also a three-leg vertex with one Higgs boson 
and two photons given by — g^Sgm„,F^5^^. 

3.3.3 MS counterterms 

We relate bare quantities to renormalized ones introducing multiplicative renormalization constants Z^ and 
(if not otherwise stated) we expand them through the renormalized S'[/(2)-coupling constant g^, 



n=l 



^' = 1 + E hfe (34) 
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where (5z|"^ are counterterms and the subscript i refers to masses, couphngs, gauge parameters and fields. 

1/2 

In particular, for all bare masses we write m = Zm rripi and for g, and cg {sg) we define p — ZpPi^, 

1 /2 

where p = g, gs, eg, sg. For a given bare field we find convenient to write (j) — (jy^i , where (jy^i is 8. 
renormalized field, and we expand through Ea. (j34p . The bare photon field represents an exception, 
and here we use 

= Z]J2 + Z]!^ Z^, = ^ ( ) , (35) 



n=l 



where and Z^ are the renormalized fields for the photon and the Z boson. Note that Z^a is expanded 
through Eq. (j34p . In addition, bare fermion fields V' (we omit flavor labels) are written by means of bare 
left-handed and right-handed chiral fields V'^ and V'"- The latter are traded for renormalized fields V'h and 
V'n expanding the renormalization constants through Eq. (l34|) . 

= i (1 ± 7^)^^, r-'- = V-^'^ (36) 

Faddeev-Popov ghost fields are not renormalized. For the bare gauge parameters introduced in Eq. ipS)) we 
use ^ = where ^ is one of the bare gauge parameters and is the associated renormalized quantity. 

Z(^ is expanded by means of Ea. ([M|) except for the case ^ — £,az, where we use 



^«-=E Tfe <l (37) 



,2 

R 

16 TT^ 
n—L 

After the expansions, we use the freedom in choosing the values for the renormalized gauge parameters and 
set ^ = 1. Next, we define a minimal MS subtraction scheme, 

<5ZW= (^-^ + A,.) AZ«, (38) 

and fix the counterterms in order to remove the poles at e = for any one-loop Green's function. The full 
list of one-loop counterterms in the MS scheme can be found in section 5 of Ref. [33j. For completeness, 
we list here all the needed results, using short-hand notations for sums over fermions [l — > charged leptons, 
u,d —f quarks), 

^/-E<., ^ = E<., ^i = E<.: (39) 

i—l i—1 i—1 

and we introduce scaled masses, Xi = mf ^/rn^ ^. In the following we drop everywhere the subscript R, 
since all quantities are renormalized ones. 

Gauge parameters. 



AZ^'^ 


1 

~ 2 


AZ«, 


Azf^ ^ 

S,AZ 






AZ^/^ 


1 

^ 2 


AZP, 


AZf^ = 


-i(A4V- 




AZ^'^ 


1 

~ 2 


AZ^^\ 


AZf^ = 




f AZ(i) 



(40) 



Gauge-boson and Higgs-Kibble fields. 



Azi^i 



23 

1 /41 1 



AZ« 



AZ'^o' 



sg_ 
3 



41 1 

~^~g 
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Masses and couplings. 



3 1 7 

icg 3 X 

2 U d 2 



1 r3 1 

[24 

23 n 



AZ(i) 



AZ 



(1) 



2X? 



X: 



x„ + + 3 (Xd + X„) 



AZ(^) 



(1) 

AA ■ 



(41) 



(42) 



Higgs-boson field and mass. 



1 

ci 



3 Xd + X, 



AZ 



(1) 



1 1 
2^ 



1 

3^' 



Xd + X, 



.(43) 



Fermion fields and masses. 



AZW=0, Azi'J=\-l + ^ 



AZW= + ^ + M 



Azi'}=-^-;- + ^ + ^-^ 



AZ, 



(1)- AZ^'^-ifl 
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1 1 
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3 5^ 
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45| 
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(44) 
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3 1 
2~4 
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QCD Counterterms. 



AZ. 



(1) _ 



2 



AZW = 



9s 2 52 



(45) 



(46) 



3.4 Finite renormalization 



We devote this section to discuss: 1) wave-function renormalization (WFR): the one- loop photon (gluon) 
and Higgs-boson WFR factors, Za (2g) and defined in Eq. fTB)) and Eq. fTT)) . are inserted in the amplitude 
of Eg. pop : 2) finite renormalization: all renormalized parameters showing up in the UV- finite amplitude 
are related to two different experimental input-parameter sets (IPSs) through the one-loop solutions of the 
SM renormalization equations. In addition, we show through a detailed analysis of finite renormalization 
that on-shell mass renormalization clashes with the simplest WST identity at hand, the doubly contracted 
relation for on-shell photons of Eq. p^ . 

In order to deal with compact expressions, we use the relations ■ pi — for on-shell photons and set 



F, 



(21). 



the amplitude for H — > 77 reads 
M 



^-^ (ei-e2)F„(0,0,-Af2) -f (ei-p2)(e2-pi)Fp(0,0,-Af2). (47) 

The expression for the H ^ gg amplitude is obtained introducing color indices and replacing g^ Sg with g^ 
and Za with Zq. 

Furthermore, after expanding the form factors Fjj and Fp and the function A4 at two loops, we split pure 
0{g^) electroweak corrections and 0((7g) QCD components, where gs is the renormalized SU(3)-coupling 
constant of the QCD Lagrangian, 

F = 



9 p(2,Ew) 

16 7r2 



9. 



16 7r2 



s ^(2,qcd) 



F ^ Fn,Fp,M. 



(48) 
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3.4.1 Wave-function renormalization 



111 the context of the MS-renormahzation prescription used in this paper, field-renormahzation constants 
are not chosen in order to compensate virtual corrections induced by WFR factors a la LSZ. Therefore, at 
variance with the on-mass-shell prescription, external legs have to be properly dressed through the formalism 
introduced in Eq.((l6l) and Eq. ^TT]) . 

Before expanding the WFR factors, we briefly recall our notation. Let S^'' be the sum of all IPI diagrams 
for a vector-boson self-energy {i — A, Z, W) or the transition between the photon and the Z boson (i = AZ) . 
We isolate tensor structures and introduce perturbative expansions according to 

n— 1 ^ ^ 

where f^" = 5^'^ — l'^" and l'^'^ = p^p'^ /p^ . The same expansion in g will be used for the Higgs boson or a 
fermion self-energy, both denoted by (i = H, f). 

The LSZ WFR factor for a photon, Z^, is fixed through Eq. ^TB]) . The transverse part of the photon 
Dyson-resummed propagator A'^ can be expressed using Eq. (105) of Ref. [12]; it involves the transverse 
parts of the photon and Z-boson self-energies and the photon/Z-boson transition. 



Ar(/) = t^-A'Ap^) + n-^AP% A'Ap^) =p'- ^4p') - ^r^^y (2y (50) 

L J p^ + m% — 2^z{P ) 

where ruz is the renormalized Z-boson mass. For a one-loop accuracy, self-energies and transitions are 
computed at order ©(g^), and factors involving the weak-mixing angle 9 are singled out in the same spirit 
of the LQ-basis formalism, thoroughly discussed in section 6 of Ref. [32] , 

si'^b') = P'sln'^l\p% = p^'-lni%^), ^ ^nW(p2). (51) 

The one-loop photon vacuum-polarization function H^^^ and the residual function T\^z^ are regular at = 0. 
In addition, because of the diagonalization procedure summarized in section I3.3.2[ also is regular at 
~ 0. Therefore, the canonical LSZ condition of Eq. fTB)) for the photon WFR factor allows to express Za 
also in the full SM by means of Ili^^ evaluated at zero-momentum transfer, as in QED, 

2- = 1 -TJ4ni'^(o)- (52) 

16 TT"^ 

Note that an analogous relation holds at the two- loop level, see section 5 of Ref. [44] and Ref. [67], where 
the same result was obtained through a background-field method analysis. 

The WFR factor for the Higgs boson, Z^, is fixed in Eq. fTT]) : respect to the photon case, here mixings 
are not present because of CP conservation. After expanding the Higgs-boson self-energy at one loop and 
introducing the real part of its derivative, E^,p, we get 

= 1 - -I^ResW (-M^), S« (P') = ^^f^- (53) 

16 TT^ Op^ 



-1 



i^Azip')]- 



Including the WFR factors of Eq. (|52)l and Eq. (|53)l in the amplitude of Eq. (|47)) , and expanding A4 through 
Eq.lO, we get finally 



A = 



^{m^'^ + [■^^''™^ + >(W (^ResW (-A^l) + 5^n«(0)) ] + ^A^(2.«co) y (54) 

Of course, an on-mass-shell prescription for the counterterms would shift sL^p and n^^-* in the one-loop 
expressions for the counterterms themselves; the associated contribution, in our notation, would be hidden 
in A^(2^'=w)^ 
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3.4.2 Finite renormalization for masses and couplings 

The second step in building the amplitude consists in performing finite renormalization; the residual depen- 
dence of the UV-finite amplitude on the renormalizcd parameters is removed through the solutions of the 
SM renormalization equations, truncated at the appropriate order. For a NLO accuracy, we need tree-level 
solutions for all 0{g^) and 0{g^ g^) terms of Eg. ([54]) and one-loop ones for the 0{g'^) component. As a 
result, the final expression of the amplitude will contain only the selected experimental IPS and will be 
ready for numerical evaluation. 

Note that the dependence of the amplitude on the renormalized parameters appears at two different 
levels: 1) explicitly, by means of the pre-factors containing the coupling constants g, sg and gs] 2) implicitly, 
through the functions M^^\ M^^.^w)^ _yy^(2,QCD)^ ^j^i^j^ 

a priori depend on the renormalized 
VF-boson, Higgs-boson and fermion masses mw, mn and m/, and the weak-mixing angle 6 (or, equivalently, 
on the renormalized Z-boson mass ruz = ruw /eg). Here the Cabibbo-Kobayashi-Maskawa matrix is identified 
with the unit matrix. 

Finite renormalization for _A4(2,qcd)^ j^W^ ^^^^ -qC^) trivially achieved by identifying the 

renormalized parameters showing up in their explicit expressions with the tree-level solutions of the renor- 
malization equations, which will depend on the chosen IPS. 

The case of M.^^\ instead, is more subtle. Let us introduce an appropriate suffix, to indicate 

that the renormalized parameters {pi,r} appear in Solving at one-loop the renormalization equations 

amounts to replacing {pi,r} — > {Pi,e + where {pi,e} are taken from experiment and belong to the chosen 
IPS, whereas {5pi\ summarize the one-loop corrections. Finite renormalization is obtained by 

Al(^) - A^(^) + Y.K1M^ ^'Z = (55) 

J J. 

For H —> 77, the dependence of Al*^^^ on g and 9 is encapsulated in the overall pre-factors of the amplitude 
of Eg. ([5^ . and the renormalized mass of the Z boson niz does not show up at one loop. Therefore, it is 
convenient to choose the on-shell masses of the W boson, Miv, the Higgs boson, Mh, and the fermions, Mf, 
to be part of any IPS, and to evaluate derivatives respect to the associated renormalized masses nin^, ttih 
and nif. 

We employ on-shell mass renormalization, as already discussed for the Higgs boson in Eq. (fT3|) . We use 
the one-loop solution of the renormalization equations for the Higgs-, W^-boson and fermion masses, 

2 2 

^ ^'^l + t|-T I^'^^- ^ (-^^b)' B = H,W, mf = Mf + ReSW (-M^), (56) 

16 TT"^ 16 7r^ ^ 

and write the amplitude of Eq. ([5^ as 

A = f^|A4(^) + ^,M^"-^^^ + f^A^«n«(0) + (57) 

16 TT^ L 16 TT^ 16 TT^ 16 TT^ J 

where we have factorized the residual dependence on the renormalized coupling constants g, 9 and gs, defining 
-^(2,Kw) ^ ^ (-^^^) + E -^iilReEW(-Af2) + 2^M/A^« ReS(i)(-M2). 

i=W,H ' / ^ 

(58) 

The last term in Eg. ([55)1 contains a dependence on g^ coming from QCD corrections to finite top-quark mass 
renormalization. At this stage, we can safely identify all renormalized parameters showing up in the explicit 
expressions for A^*^^\ A4''^'^^\ J^('^-'^'^°') and H^^-* with the tree-level solutions of the chosen renormalization 
equations. 

For electroweak corrections to _ff — > gg a similar decomposition holds. 
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provided with color indices. Note that, concerning H gg, the derivative with respect to the Higgs-boson 
mass in Eq. (j58p vanishes, since there is no dependence on ttih at LO. 

We complete finite renormalization by relating the renormalized coupling constants 9 and gs, collected 
as simple pre-factors, to three additional experimental input data. For gs, we write g| = 47ras(/^|), where 
as is the strong-coupling constant and ii^ the appropriate renormalization scale. For g and 9, instead, we 
can select two data among Gp, the Fermi-coupling constant, a, the fine-structure constant, and M^, the 
on-shell Z-boson mass. 

Finite renormalization for H —> 77.' IPS 1. The first IPS is {Gp,a\. We start writing .g^Sg — e^, where 
e is the renormalized electromagnetic-coupling constant, and use the relation for one-loop electric-charge 
renormalization. 



4 TT a 



1 



4;r 



n«(0)' 



(60) 



where T\''a \ the one-loop photon vacuum-polarization function, is defined in Eq. (|5ip . After using charge 
renormalization in the amplitude of Eq. (|57p . we observe that all terms containing YI^a' cancel out and we 
obtain 



47r I 



9 



— (2,BW) 



lTT ^. 16 TT^ 

For g, we employ the relation following from muon decay |68| 



47r 



(61) 



(\/2G,)'^'{l 



4V2 7r2 



ReE«(-M^) - S«(0) - MlSa] }, 



(62) 



where, following Ref. 



we split universal hard corrections to the muon lifetime, encapsulated in the 



W^— boson self-energy , and process-dependent components, summarized by the quantity Sq, whose 
explicit expression reads 



6 



7_44 
2sl 



Note that sg and ce are consistently fixed at the lowest order in perturbation theory by 

c2 



4 



1 



9^ 



V2GpMl ' 



(63) 



(64) 



and we are avoiding any reference to the on-shell Z-boson mass . The amplitude will finally read 



27r V J i 2V27r2 V 



4tt 



-)}, (65) 



where the total contribution from finite renormalization reads 



M 



(2,EW) 
IPSl 



ReS]W(-Ae) - l]i^^(0) 
Ml 



- Sc 



E 

i=W,H 



M[]\Re4'\-Mf) + 2j2MfM^Rei:f\-Mj). 



(66) 



Finite renormalization for H ^ 'yj : IPS 2. The second IPS is {Gp,Mr^}. We start replacing Eg. ([62]) in 
Eq. ((57|) . removing g. We get 



3/2 



.f2 r 



2V2 7r2 



+ ,2^i:i^_^(i)n(i)(o) + ^Mlm 



{2,Bw) , 3^(1) ^ReE^;^(-Af^)-Si;'(0) 



Ml 



2V2 7r2 



Q's(Mh) ^^(2,qcd) 

47r 



(67) 
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Next, we write = 1 — and use the relation between renormalized masses cg = m„/mz. After using 
vector-boson on-shell mass renormalization it follows that 



1 



Ml 



r2 ) \ 



1 



2\/27r2(M2 -M2J 



M,2 



M 



/f2 



(68) 



After employing Eq. ([68| in Eg. ([67|) . the amplitude will then be expressed through Gp, M^^, and M^; 5 a 
is given in Eg. (|63|) . sg and eg are fixed at the lowest order in perturbation theory by Sg = 1 — Cg and 
Cg = M'^/M'^, and any reference to the fine-structure constant a is avoided. 

Finite renormalization for H ^ gg. Concerning finite renormalization for H ^ gg, we start from 



A 



■M 



(2,EW) 



47r L' " ' 167r2 
and we use the result of Eg. ([62|) for muon decay to express g through Gf 



(69) 



3.4.3 The doubly contracted WST identity for _ff — > 77 at two loops 

Having introduced WFR factors and finite renormalization, we can now have a closer look at the doubly 
contracted WST identity with two on-shcU photons of Eg. (fT9|) . As previously stressed, the computation 
of the identity reguires Higgs-boson mass renormalization as a key ingredient, because bosonic couplings 
proportional to the MS-renormalized Higgs-boson mass, mn, appear at the one- loop level. Since the Higgs 
boson is emitted from a physical source, the associated momentum is on the mass shell, = —M'^, and 
rrin has to be consistently traded everywhere with M^, including radiative corrections. 

We write WTdc of Eg. (fTO|) through the same decomposition already used for the amplitude in Eg. ([TT]) , 



Wide (0,0) 



^3 „2 

Sa -^-1/2 



F^{Q,0,-Ml) - _« F,(0,0,-Af2) 



(70) 



16 7r2"" 2 

where Fp = F): ' . Next, we expand at one loop the Higgs-boson WFR factor employing Ea. ([S51) . and 
the function X through Eg.®, for = A", 

„2 r 1 -, „2 

(71) 



WX,e(0,0) = I^Ia-W + [^-(2^-") + \x^^^)Kci:%{-Ml)] + ^A'(2.Qco)l 

16 7r^ L Idtt^ L 2 f h ^ 16 tt^ J 



Since the identity holds for an on-shell Higgs boson, we perform finite renormalization employing Ea. (|55p 
and obtain 



3,2 



WX,e(0,0) = |-^ A-W+.l 



167r2 



167r2 



^(2,EW) 



ReYs^^li-Ml) -f A-;;^ ReI]^^^(~Af2 



(1) 



I 9s y(2.QCD) 

^ 167r2 



(72) 

Note that we do not need to employ other renormalization eguations, since the identity has to be proven 
both for renormalized and experimental masses and couplings. We compute the identity at one loop and we 
prove that X'^^^ = 0. The final form for the identity at two loops reads 



WX,,(0,0) = |!4|t|^ + ReE«(-M^)l + ^ A-^^.^-) | 



(73) 



Finally, it is harmless to replace renormalized parameters with the tree-level solutions of the renormalization 
equations, 



WI..(0,0) = (^/2G.)'''{^^ [a-(-™) + ReS«(-M^) 



2V27r2 

We prove that Wide 7^ 0, explicitly we obtain 



47r 



^(2,qcd) 



WXde(0,0) = (V2G^) 



3/2 



16 TT^ 



Co {-Ml , 0, ; M„ , M„ , M„ ) Im sL^^ {-Ml ) 



}. (74) 



(75) 



where Co is the scalar three-point function. The analysis of this paradox - violation of WSTI - will be 
postponed till section [ 
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4 Manipulating Feynman integrals after generation of Feynman diagrams 



In this section we summarize the techniques used after the generation of Feynman diagrams, performed 
with the FORM [39 program QraphShot [40J, the projection of the amphtude onto the form factors 
of Eg. ljlip with the projectors discussed in appendix [B] and the standard operations concerning the Dirac 
algebra. Before attempting the semi-analytical or numerical evaluation of (pseudo-)observables, we per- 
form three kinds of simplifications and symbolic manipulations: firstly, we remove reducible scalar products; 
secondly, we employ integration- by-parts (IBP) identities HT (see also Ref. [SD|) for simplifying tadpole 
diagrams; finally, we achieve an optimal level of symmetrization for loop integrals. 

4.1 Reduction of scalar products 

Given a generic loop integral, it is possible to assign a one-to-one correspondence between a specific 
reducible scalar product containing at least one integration momentum and a particular propagator. This 
simple observation allows to recursively write loop-momentum dependent scalar products in terms of propa- 
gators, and employ a basic cancellation mechanism, when non-trivial numerator structures appear, through 
the algebraic master relation 

2q-p _ 1 1 

{q'^ + ml)[{q + p)^ + m'^] q^ + m\ {q+pY+m\ {q^ + m^)[{q + pY + ni^ ' ^'"^ 

Note that the cancellation of a scalar product is associated with the disappearance of lines in the corre- 
sponding diagram; as a result, each diagram generates a set of child diagrams with a smaller number of 
propagators. 

The number of independent scalar products involving loop momenta can exceed the number of propaga- 
tors in a given diagram; therefore, some irreducible scalar products cannot be removed from the numerator 
functions. This is the well-known obstacle in achieving a full reduction for two-loop diagrams; for the special 
case of two-point functions, it can be by-passed through a judicious sub-loop reduction, as shown by the 
authors of Ref. [70] . Obviously, for a fixed number of loops and external legs, diagrams with a large number 
of internal lines, and thus more propagators, exhibit less irreducible scalar products than diagrams with a 
small number of internal lines. Note also that the choice of the scalar products which are considered as 
reducible and irreducible is to a large extent arbitrary. 

For two-loop three-point functions we have two independent external momenta and two integration 
momenta; therefore, we have to deal with 7 — 1 irreducible scalar products, where /, with 4 < / < 6, is the 
number of internal lines, and a full reducibility using Eg. (17^ is clearly not at hand. 

As an example, let us consider the F^-family diagram shown on the left-hand side of Ea. ((77)l . where 
the Higgs boson couples to photons through a couple of W bosons and a top-bottom loop. After acting on 
the diagram with the projector Pj^'^ as described in appendix |B] and saturating all free Lorentz indices, we 
remove all possible reducible scalar products and obtain: 




Here the non-planar diagram on the right-hand side corresponds to the configuration defined in Fig. [57] 
(the dot-line denotes a light fermion), where the tensor structures have been stripped and explicitly collected 
as an overall factor in square brackets; the symbol Ch is proportional to Ch pi ■ p2 + — and we 
have denoted by "reduced diagrams" all diagrams with at least one internal line less. The scalar product 
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Pi ■ qi survives as an irreducible one, showing the presence of scalar, vector and tensor l^"-type integrals, 
which appear in the combination of Ea. (|77p . 

The reduction procedure is repeated iteratively on the reduced diagrams; as a matter of fact, the V" 
topology generates a set of eight sub-topologies, including simple factorized topologies (products of one-loop 
functions) and two-loop vertices, self-energies and vacuum diagrams, as illustrated in Fig. [6l 

During the reduction procedure, we devote special care to preserve the canonical routing of loop momenta 
defined in appendix [Al In Fig.[7]we show an explicit example for reducing the scalar products of a y^-family 
diagram. Note that this configuration can appear after performing the reduction of the V" diagram using 
Eg. ([77)1 . In step (la), the scalar product q2 • pi is removed; this operation leads to the appearance of an 
additional integral belonging to the family which does not possess the standard routing of momenta as 
defined in Fig. The canonical routing of momenta is recovered in step (lb) through a shift in the loop 
momenta which acts also on the remaining scalar product qi ■p2 and generates a new reducible scalar product 
12 ■ P2- The latter is then removed in step (2) , leading to an additional integral of the family. 

In some cases, only the scalar configuration survives after the projection procedure. An example is the 
configuration in the left-hand side of Eq. ((78| , where the Higgs boson couples to two photons through a 
couple of Z bosons and a top-quark loop. After applying the projector P^'^ and removing reducible scalar 
products, one obtains the decomposition 




A A 



Figure 6: Generic child topologies of the V" parent topology. The five-line diagram is obtained by removing one 
line of the diagram; the second line contains the child topologies of {V^ , S"^ and B x B). The third line 
contains the topologies , B x A and T^, obtained by removing one line from the diagrams above. The arrows 
indicate the correspondences between parent and child topologies. 
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Pa 




Figure 7: Reduction of scalar products for an example of the family. The equations are valid for Pi = pi = 0. 




where the symbol Ck is proportional to Ck cx 32 (w^ + v'i) {pi ■ p2 - M| + 2 Mf) - 128 v+v- {pi ■ 
P2 + 2Mf). Here, v± are the V ± A couplings Ztt. Clearly, here no irreducible scalar products for the 
topology remain, and only the scalar configuration survives. 



4.2 Vacuum diagrams 

The removal of scalar products illustrated in Fig. [6] shows that the vacuum diagram appears after the 
reduction of the S'^' self-energy. The explicit evaluation of T'^ has been carried out in Ref. [71] . In addition, 
if the external momentum of the S'^ integral corresponds to one of the photon momenta, the S'^ integral is 
a vacuum integral with an increased power of one propagator. Furthermore, if additional loop-momentum 
dependent scalar products appear in the numerator of the S'^ diagram, the integral has to be expanded in 
powers of the external momentum around zero; therefore, higher increased powers of the propagators can 
arise. 

Vacuum integrals with increased powers of the propagators can be related to the tadpole integral 
using the traditional IBP method. In particular, the relation between , the tadpole diagram with one 
increased power of the third propagator (shown diagrammatically through a dot on the corresponding line) 
and is given in Eq. ((791) by 

(2 - n) (QO) 




(n-3) 



A(to^, 7712, mg 
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fY^ + {n-2) 'jr^'i^rK^ G^, (79) 



with \{x, y, z) — + + — 2xy — 2xz — 2yz. 

After projecting the amphtude on the relevant form factors, it turns out that only the tadpole diagrams 
and , related to each other through Eq. ([79|) . contribute to the amplitudes ^ 77 and H ^ gg. 



4.3 Symmetrization of loop integrals 

It is essential to exploit the symmetries of each diagram in order to reduce the number of integrals to be 
calculated and to identify equal configurations. We have taken into account the symmetries of the appearing 
one-loop and two-loop topologies, summarizing them in Tab. [Hand Tab. [2l In both tables, the first column 
denotes the topology and the second column enumerates the different symmetry transformations for a given 
topology; the third column contains the transformation of the loop momenta qi and q2 and the fourth column 
the corresponding interchange of masses and external momenta. The identity transformation and the total 
reflection of all external momenta {pi — > — pi, P2 ~* ^P2 and P —P, corresponding to the loop- momentum 
transformation qi — > —qi and q2 —^72)) which leaves the loop integral unchanged, are not listed in Tab.[l] 
and Tab. [H The largest number of symmetries can be observed for the family. 
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q - 
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— > m2, m2 - 




ms 


^ mi. 












Pi - 


^ P2, P2 


-P, - 


P - 


>Pi 



Table 1: Symmetry transformations for one-loop topologies. 



As an explicit simple example concerning the application of symmetry transformations, let us consider the 
case of the 1^^-family integrals arising after the reduction of the V^- and V^^'-type Feynman diagrams shown 
in Fig. [8] (a) and (c). Note that the momenta routings of the two child y ^-family diagrams in Fig. [8](b) and 
(d) are different; the application of symmetries aims to map the various integral representations onto a single 
one in order to allow for cancellations. The family has three basic symmetry transformations, summarized 
in Tab. [2] and graphically illustrated in Fig. [9l here the first diagram corresponds to the standard integral 
representation of V as defined in appendix El the symmetry transformation (I) amounts to exchanging the 
first and the second line of the self-energy insertion in the diagram; (II) corresponds to an exchange of the 
fourth and fifth lines (in this case also the external momenta are interchanged, pi —P)', symmetry (III) 
is a combination of symmetries (I) and (II). The two V child topologies in Fig. [Sfb) and (d) are related 
through the symmetry transformation (III). 

Finally, we stress that the procedures of reducing scalar products and symmetrizing loop integrals have 
been recursively performed in order to achieve a maximal simplification of the amplitude. Next, the remaining 
integrals are classified as scalar-, vector- and tensor- type integrals according to the number of irreducible 
scalar products in the numerators. The Lorentz structure of the loop integrals are expressed through the 
external momenta pj, with j = 1, 2, and the metric tensor introducing suitable form factors. A FORTRAN 
code is generated and the form factors are then evaluated numerically employing the NAG library [72] , 
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Table 2: Symmetry transformations for two-loop topologies. 



5 Behavior of two-loop diagrams eiround a normal threshold 



In this section we will discuss one of the main results of the paper: by carefully studying all singularities 
of Feynman diagrams we propose the two-loop implementation of a renormalization scheme which cures 
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anomalous behaviors of the amphtude. 

Feynman diagrams have a comphcated analytical structure as functions of the external Mandelstam 
invariants and the internal masses. A frequently encountered singular behavior is associated with the so- 
called normal thresholds: the leading Landau singularities (73) of self-energy-like diagrams which can appear, 
in more complicated diagrams, as sub-leading singularities. In this section we discuss how the amplitudes for 
H 77 and H gg behave around a normal threshold, with special emphasis to the problem of possible 
square-root and logarithmic singularities. 

Without loss of generality, let us consider the case of the ^ 77 decay in the setup where all light 
fermions are taken to be massless. As far as normal thresholds are concerned, we have the possibilities 
= , , 2M^ , 2M^ , 2Mt, as illustrated in terms of cut diagrams in Fig. [TOl Note that, as 
observed by the authors of Ref. [28], there is no cut at M^^ = even in presence of massless fermions; indeed, 
the two-particle cut of the first diagram of Fig. [10] is zero because of the helicity structure of the diagram. 

The = M„, , cuts are identified by the configuration shown in the second and third diagrams of 
Fig. [To] Note that, for H 77, the imaginary part is not exactly zero below the single- ly threshold; this 
is due to the introduction of complex masses for vector bosons, as we will explain in detail in section 15.31 
It is worth mentioning the behavior of the one-loop amplitude which is rapidly decreasing for small values 
of , as shown in Fig. [Tl] for H gg. Here, once again, we use a complex W mass and, therefore, the 
imaginary part is different from zero even below the it threshold. We have performed a dedicated analysis of 
the behavior of the amplitude around single thresholds which shows that regular behavior is a consequence 
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Figure 10: Sample two- and three-particle cut diagrams for H 77. The mass of the up and down quarks is 
neglected. 



of a delicate mechanism of numerical cancellation among several diagrams, e.g. V" , doubly coUinear and 
simply collinear. 
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Figure 11: Real and imaginary parts of the one- loop H ~* gg amplitude. 

Having discussed the single vector-boson thresholds, we move to the double ones, identified for example 
by the cut of the fourth diagram of Fig. [TOl 

5.1 Square- root singularities 

In this section we are interested in the problem of possible square-root singularities of the amplitude. 
When present, they are unphysical, although integrable; the solution to the apparent puzzle consists, as we 
will argue, in replacing real masses of unstable particles with their complex poles. 

In order to illustrate the problem, let us consider the following integral, related to a generalized one-loop 
two-point function: 



dx X 



X + fi — iO 



(80) 



where n is a non-negative integer. When (3 = 0, with = 1 — the integration contour is pinched 
between the two singular points a;± = (1 ± /3)/2, and the integral of Eq. ((80l) is singular, with a branch point 
of the two-particle cut = 1/4; for a > 1 we find 



2/ \2' 2 



(81) 



where B is the Euler beta function. To be more explicit, let us consider the UV decomposition of section [2] 
for the generalized one-loop two-point function with equal masses 



(82) 
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where k and / are the powers of the two propagators. Concerning the finite parts for the fc = ^ = 1 and 
k ^ 2,1 ^ 1 configurations, we obtain 

ui^ 1 
Bo(l, 1 ; s,to,to; 0) = 2 - hi — — - - (3Lj3, Bq{2, I ; s,m,m - 0) ^ Lp, 

s + tu p 

In general, for fc + ? > 2, we have 

Bo(fc,^; s,m,m; 0) /? ^ 0, (84) 

and we can easily conclude that the generalized one-loop two-point function has an (unphysical) singularity 
at s = Am? yk,l \ k + I > 2. For practical applications, one has to distinguish between the above- and 
below-threshold regions, 

— In i+« 



/^^ = -a^ < 0, L/3 = i (arctan-r^ + n] ^ in, s ^ Am^ \^ . ^ ' 



^a^ >0, Lf3 = lnj^ - in in, s 

2a 
1-a 

Therefore, the generalized one-loop two point function with k + 1 = 3 generates a square-root l//3-divergent 
behavior, with = 1 — 4M^/M^ or = 1 — 4M^/M^, associated with the normal threshold related to 
the two-particle cut (the leading Landau singularity). 





contraction 



Figure 12: One-loop vertex with two equal masses and its sub-diagram giving the sub-leading square-root singularity. 

Let us now consider three-point functions, and investigate the behavior of the vertex of Fig. 1121 around 
the normal threshold located at s = Am?. The leading Landau singularity of the vertex is the so-called 
anomalous threshold; the normal threshold shows up as a sub-leading singularity. Since the sub-leading 
singularity for a graph is the leading one for any of the contracted sub-graphs, we easily conclude that the 
singular behavior of the one- loop vertex of Fig. [T^] around s ~ A m^ is due to the sub-graph where we shrink 
a line to a point; the singularity is a branch point in the complex s -plane. 

The same argument can be repeated for all diagrams with any number of external legs where we can cut 
two and only two lines with mass m; any normal threshold will be a sub-(sub- . . . ) leading singularity, and 
a 1//3 behavior will show up only if the reduced sub-graph, responsible for the singularity, can be reduced 
to the generalized one-loop two-point function of Eg. ([55]) . Bq{s, to, to) = Bq(2, 1 ; s, m, m; 0). A Bo{s, to, to) 
function is related to the derivative of a Bo{l, 1 ; s, to, to) function and therefore emerges in the computation 
whenever we include WFR factors for the external legs or we perform finite renormalization of the mass of an 
internal particle, as depicted in Fig. 1131 In the second case, generalized one-loop triangle functions appear, 
and it is known that they can always be reduced to Bq functions using IBP identities (this is just another 
way to say that the normal threshold is a sub- leading singularity for a generalized three-point function) . 

Concerning genuine two-loop diagrams, we observe that a Bq configuration can only arise if we have a 
self-energy insertion in a two-loop diagram. A two-loop vertex containing a self-energy insertion, leading 
to a l//3-divergent behavior, is depicted in Fig. 1141 For this diagram it is possible to find a representation 
where the singular part is completely written in terms of one-loop diagrams, as shown in the figure. The 
remainder can be cast in a form suited for numerical integration. 

Note that the unphysical l/f3 behavior, generated by the diagram of Fig. [TU exactly cancels the one 
coming from mass finite renormalization. This is strictly true only in a complex renormalization scheme 
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Figure 13: Singular (3 ^ behavior at the normal m threshold coming from WFR and mass finite renormalization. 
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Figure 14: Singular 13 ^ behavior at the normal in threshold coming from self-energy insertions {V" topology). 
The generalized one-loop triangle function in the right-hand side is related to the generalized one-loop self-energy of 
Eq.((83l) through IBP identities. 



because, for on-shell real masses, only the real part of the self-energy in Fig. [13] is taken and the cancellation 
does not take place for the corresponding imaginary part generated by the two-loop diagram of Fig. 1141 
This is not surprising at all: self-energy insertions, signaling the presence of an unstable particle, should 
not be there. They are the consequence of a misleading organization of the perturbative expansion; Dyson- 
resummed propagators should be used and complex poles should replace real on-shell masses. The remaining 
1/(3 singularity is therefore coming only from the wave- function renormalization of the Higgs boson. This 
unphysical behavior is strictly connected to the problem of defining a proper WFR for an unstable particle, 
as pointed out in section [3. II 

It is interesting to note that for the ti threshold in H gg the Bq -functions that are potentially 
dangerous always appear multiplied by as it happens for QCD corrections; the same is not true for 
pseudo-scalar Higgs decay, cfr. Fig. 4 of [9]. 

Finally, we observe that a more severe behavior associated with (3^0 should not show up; for instance, 
we have verified that l//3^ terms which appear as a consequence of the reduction procedure for the iJ — > 77 
decay are of the form F{Ml)/{Ml - 4M^), with i^(4M^) = 0. 

5.2 Logarithmic singularities 

Let us consider the scalar two-loop diagram of the family shown in Fig. [151 ^^nd derive the corre- 
sponding integral in parametric space. We introduce the quadratic forms 

/ 1 \ ^ 1 1 4 

X(x)=(^x--j --/3^ ax,y) = x{x^l)y^ + -{l-P^), p^^i-^^ (86) 

where s = —P^ > and m is the mass of the solid line (wavy lines correspond to massless particles). We 
obtain 

Jo y x{x) 



Li2 (1 - ) _ Lis (1 - ) 
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Figure 15: The irreducible scalar two-loop vertex diagram of the 1/^ family showing up a logarithmic divergence. 
Solid lines represent a massive particle with mass m, whereas wavy lines correspond to massless particles. 



Since we are interested in the behavior around /? — > 0, we spUt into a singular and regular part, 



VJL = ^ 



dx dy 



The singular part V^f will be written as 74j 
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where we have introduced the short-hand notations 
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with T = {\ — t)y + t and T = t/{l — t) > 0. /(/:) can be split into two parts, 
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(89) 
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(91) 



Here Xi = —X, X2 = 1 — X and B{x,y) is the Euler beta function. The second term of Eq. (PT|) is regular 
for /? = 0; the first term, instead, shows a singularity due to the fact that A ~ (1 — y)^ for f3 ^ 0. However, 
we have a singular behavior only if < X < 1, which requires y > ymin = max{0 , {t — l/2)/(t — 1)}. Being 
interested in the leading behavior for /? ^ 0, we can extend the integration domain in the first term to [0, 1] 
without modifying the divergent behavior of the diagram. The singular part is then given by 
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where < Res < 1/2 and Fi denotes the first Appell function. In order to obtain the expansion corresponding 
to /3 — > 0, we close the integration contour over the right-hand complex half-plane at infinity. The leading 
(double) pole is at s = 1/2; therefore, we obtain 

J(i) = -i ln(-/32-iO)-hC'(l), /3^0. (93) 
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Inserting the result into Eq. ((89| and using dtr^/^ (1 - t) hit = - 27r, we get 

^sf„g = ^ln(-/?'-*0)+O(l), (3^0. (94) 

If the massive loop in Fig. [15] is made of top quarks, the contribution of the integral to the amplitude 
behaves like /S^ and, therefore, the logarithmic singularity is /J^ -protected at threshold; however, the 
same is not true for a VF-loop. Our result of Eq. (|M|) is confirmed by the evaluation of of Ref. [TD] in 
terms of generalized log-sine functions. Starting from Eq. (6.34) of Ref. [10] and using the results of Ref. [75] 
we expand around 9 = tt, where x = e'^ = (/3 — 1) (/3 -t- 1), with < 6 < tt. This gives for the leading 
behavior of below threshold ln{d — tt), where \n{—(3'^) = ln{9 — Tr)^ — In 2. The same behavior can 

also be extracted from the results of Ref. [9] . 

Logarithmic singularities of the kind discussed in this section are a remnant of the one-loop Coulomb 
singularity of one-loop sub-diagrams. An alternative approach that automatically resums large Coulomb 
singularities at threshold has been pursued in Ref. [76]. The reader should be aware that in the pseudo- 
scalar decay the real and imaginary parts of the form factor may be significantly different from the lowest 
order perturbative one|l|- 



5.3 Complex masses 



In this section we set up and discuss our implementation of a consistent and gauge-invariant treatment 
of unstable particles in NNLO radiative corrections. 

Our two-loop renormalization scheme has been described in details in section [3] where counterterms have 
been introduced, different choices of IPSs considered and finite renormalization of lagrangian parameters 
discussed. In short, this represents the so-called 

• RM - scheme 

where masses are the real on-shell ones; it gives the extension of the generalized minimal subtraction scheme 
up to two loop level. The analysis of section [5] has shown the presence of pathological features and the 
cure proceeds in two steps. Our first, pragmatical, solution to the problems induced by unstable internal 
particles has been presented in Ref. [38j (for an alternative approach to the problem of unphysical threshold 
singularities connected to WFR factors, see Ref. [77]); the corresponding scheme will be termed minimal 
complex mass scheme (hereafter MCM), the first emergency kit. 

• MCM - scheme 

To evaluate the amplitude we start by removing the Re label in those terms that, coming from finite 
renormalization, violate WSTIs. For instance, when we compute the doubly contracted WSTI for the full 
two-loop amplitude in ^ 77 we obtain the result of Ea. (|75|) : pure two- loop contribution to the WSTI 
gives T,^h\—M'^) while finite renormalization gives its real part KeI^^H\—M'^). Therefore, the WSTI is 
violated above the WW threshold, as shown in section [3.4.31 

Furthermore, we decompose the amplitude for _ff ^ 77 according to 
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(95) 



(2) (2) (2) 

and prove that, as expected, A}^' ^ A)^ and Al-cm separately satisfy the WSTI. The latter fact allows ns to 
minimally modify j4„7, by working in the complex-mass scheme of Ref. |60j : we include complex masses in 
the gauge-invariant leading part of the two- loop amplitude as well as in the one-loop part. 

The decomposition of Eg. ([95)1 deserves a further comment. As we stressed in section [5TT1 there are three 
sources of terms: a) pure two-loop diagrams of the V^, family, i.e. bubble insertions on the internal 
lines of the one-loop triangle; b) W-uiass renormalization, i.e. on-shell W self-energy x the mass-squared 



discussion with M. Spira on the last two points is gratefully acknowledged. 
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derivative of the one- loop W triangle (the latter giving rise to c) Higgs wave- function renormalization 

factor X lowest order (the former giving rise to 1//3). One can easily prove that only c) survives in MCM 
and a,b), which are separately singular, add up to a finite contribution (/? 0); their divergence is an 
artifact of expanding Dyson-resummed propagators in the on-shell approach. 

The In /3 -dependent term originates from pure two- loop diagrams of the family and it is a remnant 
of the one-loop Coulomb singularity of one-loop sub-diagrams. 

The amplitude for H gg is different in some points, for instance we have no diagram with a 
logarithmic behavior. Also the violation of WSTIs, described above, is specific to _ff — > 77 since we have no 
one-loop bosonic triangle. Further details have been presented in Ref. [75] . 

The one-loop ^ 77 amplitude, with a complex W mass, is shown in Fig. [TBI around the WW threshold 
including a comparison with the real- -mass amplitude. We have also analyzed the effect of (artificially) 




145 150 155 160 165 170 145 150 155 160 165 170 

M„[GeV] M„[GeVl 

Figure 16: Real and imaginary parts of the one-loop ff ^ 77 amplitude with real and complex -boson mass. Note 
the sizable difference with Fig. 1111 

varying the imaginary part of the W -boson complex mass (more details can be found in Ref. ^78] ) , showing 
that our complex result reproduces the real one (with a complex Higgs-mass renormalization condition) in 
the limit F,,,- — > 0. Here, a comment is needed: we introduce complex masses as poles on the second Riemann 
sheet, a fact that requires a careful analytical continuation of loop integrals. 

Let us consider the decay ^ 77 in the RM scheme; on-shell renormalization spoils the complete can- 
cellation between pure two-loop diagrams and H^-mass renormalization mentioned before. As a consequence 
the RM scheme badly fails to approximate the complex mass results above the WW threshold; here RM is 
not the smooth limit F„. — > of MCM and we understand the reason, RM is missing a cancellation which 
is instrumental in building a consistent theory of unstable particles. If RM scheme is further modified by 
arbitrarily forcing this cancellation the corresponding result is the smooth limit just mentioned. Thus we 
have an additional argument to reject RM; at two-loop this scheme is wholly inconsistent. 

It is worth noting that cancellation between mass renormalization for some internal line and bubble inser- 
tion in the same line, as far as divergent terms are concerned, is a strict consequence of Dyson resummation 
with complex poles, see Eq. (195) of Ref. [H] and consequent discussion. 

Finally, we mention the fact that MCM can be forced to agree with RM in the limit F^i, only if we 
adopt the unjustified approach of continuing the W self-energy in different Riemann sheets depending on its 
origin, irreducible two-loop diagrams or finite renormalization. 

In a nutshell, the MCM scheme has been designed to cure the unphysical infinities of two-loop amplitudes, 
namely those points where the amplitude is artificially infinite; it does not deal with cusps associated with 
the crossing of normal-thresholds present in Arcm, as described in detail in Ref. [78]. Note that this is not 
only an aesthetical issue but also a concrete problem in assessing the impact of electroweak NLO corrections 
on, say, Higgs production via gluon gluon fusion: observing the effect of NNLO QCD corrections with 
respect to NLO ones one is lead to understand possible sources of additional large corrections. Electroweak 
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corrections, typically around the WW threshold, can reach a 10% in the MCM scheme due to a magnified 
cusped behavior. 

• CM - scheme 

Further to the last point, we have undertaken the task of introducing the (complete) complex- mass scheme 
(CM), based on Ref. [60j, as explained for a two-loop calculation in Sect. 10 of [44]. This means that all 
two-loop diagrams must be computed with complex , masses, whereas the top quark mass is the 
on-shell mass; note, however, that we keep the external Higgs boson on-shell and do not perform the ultimate 
step of introducing a (complex) pole residue and the associated partial width for the decaying Higgs boson. 

For _ff — > 77 (in gauge) we have a bosonic triangle at one loop, which contains a factor uih (from the 
H — (f> — (f) vertex, see Fig. |4]), which needs to be renormalized through the replacement of the renormalized 
Higgs mass with the physical Higgs-boson mass. This fact introduces the real parts of Bq functions, which 
lead to a violation of the WST identities above the WW threshold, as described in Eq. ((75|) (note that 
wave- function renormalization factors never pose similar problem); also in the CM scheme the corresponding 
real label is removed, even if the external Higgs boson is assumed to be an on-shell particle (see comment 
above) . 

As a final note one can say that CM scheme is the default for our results, the other schemes being assigned 
to the role of benchmark. 



6 Extraction of collinear singularities 

In this section we discuss the problem of collinear singularities showing up in the calculation of a given 
(pseudo-)observable. Any method that aims to produce theoretical predictions for (pseudo-)observable quan- 
tities organizes the calculation of the corresponding ^-matrix element into several building blocks, and the 
analytical structure of the total amplitude will not necessarily be the same of the single components. From 
this point of view, collinear singularities are a clear example: sometimes, a collinear-free amplitude is split 
into components which are separately divergent in the collinear regime. Therefore, singularities must be 
regularized and singular terms have to be extracted. 

It is worth noting that no numerical evaluation can be attempted before two basic steps have been 
performed. On the one hand, all singular terms (ultraviolet, infrared, collinear) of the amplitude have to 
be extracted, and their cancellation or absorption into parton distribution functions have to be explicitly 
checked. On the other hand, all enhanced terms have to be isolated, such that numerical integration is only 
limited to smooth remainders; at this level there is no need to worry about the length of the remainders. 

Focusing our attention on electroweak processes, several methods aimed to deal with collinear singularities 
have been developed and presented in the literature, as in Ref. [79j . for one- loop leading logarithms in 
electroweak radiative corrections, and in Ref. [SOji for two-loop electroweak NLO logarithmic corrections to 
massless and massive fermionic processes. 

Concerning electroweak corrections to the decay of the Higgs boson into two photons or two gluons, 
collinear divergencies are related to the coupling between photons or gluons with light fermions. Therefore, 
there is no substantial difference between the two processes and in the following, without loss of generality, 
we will concentrate on the process H 77. On the contrary, the QCD corrections generate special types of 
divergencies in the gluonic decay, because of the three-gluon coupling. 

A common approach to the problem of collinear divergencies is to consider all light fermions of the 
theory as massless states; dimensional regularization is then used to control the collinear behavior of single 
components of the amplitude. In our approach, we prefer to keep the physical light-fermion masses to act as 
regulators, and to express the collinear behavior in terms of logarithms of these masses. After analytically 
checking that singular parts cancel in the total, we can safely get rid of the regularization parameter and 
include all collinear-free remainders into the total amplitude. These finite parts will be cast in a form which 
is functional to numerical integration. 

At two-loop level we encounter three different situations illustrated in Fig. [iTl 1) one and only one of 
the two external photons is coupled to a light-fermion current (second, third diagram and fourth heavy-/' 
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diagram); 2) both photons are coupled to the same light- fermion current in a loop (first diagram); 3) one 
photon is coupled to a light-fermion current in one loop, the other photon in the other loop (fourth diagram 
with light /'). 



7 7 7 7 




Figure 17: Complete list of diagrams with photons coupled to light fermions (/). In the figure /' denotes any fermion 
(light or not, equal to / or not). The capital B indicates a boson which can be either & Z or a W . 



The treatment of the above-mentioned configurations simplifies when taking into account the reduction 
(g) symmetrization procedure described in section 31 This aspect is considered in section [6.11 where we show 
that all diagrams of type 2) cancel. Diagrams of type 1) and 3) are then discussed in section and in 
section 16. 3[ where we analytically extract the coefficients of the coUinear logarithms. In section 16.41 the 
special case of the V^' configuration is considered; divergencies which are peculiar of the QCD corrections 
to H ^ gg are shortly reviewed in appendix ICl 

6.1 Collinear behavior and tensor reduction 

In this section we address a general question which has important consequences in classifying those 
collinear configurations which are actually needed for the calculation. For the ^ 77 process, the diagrams 
where photons couple to light fermions belong to the and V" families with the mass patterns depicted 
in Fig. 1171 Two naive approaches to the treatment of collinear singularities would consist on the one hand in 
evaluating all the possible tensor integrals associated to these families; on the other hand, to express them 
in terms of some set of Master Integrals (hereafter Mis) using IBP identities. 

The first option is not very convenient, since it requires to evaluate explicitly a large set of integrals. 
However, also a blind application of reduction through IBP identities has a drawback: the Mis representing 
the basis integrals cannot be a priori predicted (e.g. by using the standard Laporta algorithm [81]; construc- 
tive approaches can be found in Ref. [52]), and complicated coUinear-divergent Mis could show up in the final 
answer. In this section we show that our approach, based on the reduction (g) symmetrization procedure of 
section m represents the optimal solution: it identifies the smallest set of collinear-divergent integrals whose 
structure is simple enough to allow for an analytical extraction of all collinear singularities. 

To introduce our argument, we start considering an -point one-loop function with external momenta 
Pi, ... ,Pjv and with scalar products in the numerator. After introducing = and = Pi + • • • +pf , 
we consider the integral 

Sn;Af)^A I ^"g {!l'-w" . H = (g + P)'+m?, (96) 

and we perform a standard-reduction procedure to simpler functions; taking for instance four-point functions, 
it is a well-known fact that 

Sn;iU)^Y. ^» + E Go{Pi,Pf) + '^yfc Do{p!,PIpI) -f P, (97) 

i i,j 

where Bq^ Cq and Dq are scalar two-, three- and four-point functions and R is the so-called rational term. 
Let us consider, in particular, the following example: 

3 

Sn;i{q-Pi) = ^ = - VDhPh, (98) 
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where the matrix H is given by Hij = — pi ■ pj, G — det H is the Gram determinant associated with the 
four-point function and Du are standard form factors [51j . In standard reduction, one goes on expressing 
the Dii form factors in terms of Dq and of three-point functions, with inverse powers of G. However, a more 
careful apphcation of the method will make use of 

Du = -lH-/d„ d, = D^^+'^^D^;'^~2K,Do, ^ ^ {P^ - P^, + m'^ - mU), (99) 

where Dq "^ is the scalar triangle obtained by removing the propagator i from the box. Therefore, we obtain 

1 ^ 1 
^n;4((7-Pi)= 2 E H-/Hudj^-di, (100) 

without explicit factors involving G. Furthermore, from Eq. ([Ml) . we see that the coefficient of Dq in the 
reduction is (toq — —p\)/2. Note that at the leading Landau singularity of the box, corresponding to the 
anomalous threshold [83], we must have 

q^+ml^O, {q + Pif+mj^O ^ q ■ p^ ^]^{ml - m\ - p^), (101) 

which is equal to the coefficient of the I?o function. This is a general property: a careful application of 
standard reduction to an iV -point function with any scalar product gives as coefficient for the scalar Ap- 
point integral the value of the scalar product at the anomalous threshold. 

To summarize, in standard reduction for a -point function each reducible scalar product in the numer- 
ator is replaced by a difference of propagators plus a K factor, predicted by factorization properties at the 
anomalous threshold. The procedure can be continued and one finds {N — 1)— point functions with reducible 
and also irreducible scalar products; for the latter inverse powers of Gram determinants will remain. 

Imagine now that our iV -point one-loop function is a sub-diagram (with loop momentum (72) of a two- loop 
diagram (with momenta gi, 92)- The numerator will contain, in general, reducible and irreducible scalar 
products. If only reducible scalar products are present and if, after algebraic reduction iV — 1 (as we 
said earlier, no inverse Gram determinants), the coefficients of the corresponding scalar, vector or tensor 
one-loop diagrams turn out to be zero, then the two-loop diagram will not appear in the final result and only 
its reduced child diagrams will do. In particular, if the original two-loop diagram is coUinear divergent, the 
singular behavior can be read off from its sub-diagrams, which is a simpler problem because one propagator 
less is involved. 

This is what happens with the first diagram in Fig. [171 as described in appendix [Bl the physical content 
of the diagrams for the process iJ — > 77 can be extracted contracting the tensor in the amplitude by means of 
the projector P^" = D^i-^^ of Eq. (l236p . After the standard reduction 61egs 51egs is apphed (see section[3]), 
when m/ — > and for arbitrary and my , we obtain 0: 




■^In this section we use low-case letters for denoting the masses of light particles, which will be neglected after proving the 
cancellation of coUinear logarithms. 
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where the coefScient is oc (n — 4) P^. Since the scalar and vector are UV finite, we can take the Umit 
n ^ 4; therefore, as we anticipated, six-propagator terms disappear from the projected and only reduced 
diagrams with at most one photon coupled to light-fermion lines survive. Concerning reduced diagrams, the 
collinear-divergent ones are always of the type and V^, collected in Fig. [TH] Note that for the type 
only the scalar configuration survives, while for V'~' we have to deal with tensor integrals up to rank two. 




Figure 18: Definition of all Mis with one external massless particle coupled to a light particle for the process H 77. 
The dot-lines indicate light particles of mass m, the wavy line a massless particle and the solid lines whatever particles. 
The mass M is strictly heavy, while Ma and Mb can also be light. 

The same procedure will then be applied to the second and third diagram of Fig. [iTl In these cases the 
coefficient of does not vanish; in the end, however, only scalar configurations survive. 




(103) 



The coeflicients and are independent from the loop momenta. The collinear-divergent reduced 
diagrams belong to the families VJ^ and VJ^ depicted in Fig. [THl the new object to be computed is the scalar 
diagram Vj,^ , also defined in Fig. [TBI 

For the last diagram of Fig.[T7]we have to distinguish two cases, with /' heavy (i.e. /' is the top quark) 
or light. In both cases, the coefficients of the tensor V do not vanish, but some universal structure in qi 
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and q2 can be identified, irrespective of the coupling between the boson B and the fermions / and /'. In the 
case of hght /', the reduced diagrams show also a doubly coUinear behavior of the type, characterized 
by a universal structure in the loop momenta. Indeed we obtain 




where the coefficients C^^, C^^ and C^"^ depend on the coupling between the boson B and the fermions. The 
diagram VJ^ is defined in Fig. [TSl while V^" and V^^ are shown Fig. [111 the reduced diagrams (r.d.) are again 
of the form VJ^ and V,"^ of Fig. M 




PI PI 



Figure 19: Definition of all Mis with two external massless particles coupled to a light particle for the process 
H 77. The dot- lines indicate light particles of mass m or m' (m' can be equal to m or not), the wavy line a 
massless particle and the solid lines a heavy particle of mass M. 



Summarizing, we can say that the reduction ® symmetrization procedure allows us to identify the smallest 
sub-set of all diagrams with coUincar divergencies, collected in Fig. [TH] and Fig. (THl they can be taken as 
Mis for the set of simply and doubly coUinear configurations. The extraction of the corresponding coUinear 
logarithms will be treated in the next two sections. 



6.2 Vertices with one photon coupled to light fermions 

In this section we extract the coefficients of the coUinear logarithms for the Mis shown in Fig. [TH] and 
compute the associated coUinear-finite parts. For the configuration (first diagram in Fig.[TH]), we use the 
result derived in section 5.1 of Ref. [47' , and evaluate the limits — > and = rii/^ = m ^ 0, 
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where Xe ~ P'^ x {\ — x) {\ — z) + p^x z [1 — x) ^- (1 — a;) -I- x and the two-loop UV factors F^^, 
with i = 1,2, have been defined in Eq.Q. We readily identify Xe with the polynomial associated with the 
Feynman-parameter representation of the one-loop two-point function Bo{l, 1, (1— z)P^ + zp2, Ma, Mb) {x is 
the Feynman parameter). In particular, the coefficient of the coUinear logarithm can be written through the 
one-fold integral representation of a one-loop function. 



V = 




l--A„.(.)--ln- 



dz 




+ Kla^, (106) 
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where A^v can be read in Eq.(IT|) and Kj^g^ denotes the colhnear-free remainder (see Ea. ljlOSp ). This simple 
result shows a feature that we will encounter also in more-complicated configurations: coUinear singularities 
can be represented through objects with well-known analytical properties, and the cancellation of all coUinear 
logarithms at the amplitude level can be analytically verified. 

It will be shown in the following that the coUinear behavior of the remaining three vertices of Fig. [18] 
is also embedded in a one-loop integration. The key observation is that all these configurations contain the 
product of two propagators of the same type: 



(g^ -|- m'^)[{q + p)^ -|- m? 



with 




9 = 91 

q = qi 

9 = 92 
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In particular, let us consider the scalar configurations V^[l] and VJ^ (scalar configurations of the second and 
third diagrams in Fig. 118^ : we define J„ as the scalar sub- loop containing g, where N denotes the number of 
g -dependent propagators in addition to those of Eq. (|107p . For VJ^[1] and Vj,^ we have N — 1 and consider 



Ji = 



d"q 



1 



(g2 + m2)[(q+p)2 +m2][(g-g2)2 + Af< 



Next, we introduce Feynman parameters z and y, integrate over q and set n — A, obtaining 

Ji^ J^dz J^dy^, V =[A-y{q2+pf]y + m^l-y), A={q2+pzf + MI 



(108) 



(109) 



The result shows that the singularity for m = is generated, in parametric space, at the point y — Q\ 
we introduce = Ay -\- m?, a simple polynomial having the same coUinear properties of V, and add and 
subtract 1/Vb at the integrand level, getting 



Ji = 



dz Az — m? 

— In 

A s 



dy 



dz 

^0 y 



A-y{q2+py 



+ O (m2) 



(110) 



where the first term is the coUinear-divergent part of Ji. The complete expressions for V£[l] and VJ^, and 
the related coefficients of the coUinear logarithms, are obtained inserting the result for Ji in the q2 integrals. 



V;?[l] = ^ln^ f'dz f 
«7r2 s Jq J 



1 m 
^^In — 





2 j-1 

dz 



[(92 + Pi ZY + M2] [{q2 + pi)2 + Af2] [{q2 + P)2 -f M2] 



dS2 



/o J [{q2 + Pzf + Ml]{ql + M^)[{q2+pif + M^][{q2 + Pf + Ap- 
with the following diagrammatical correspondences. 



+ Kc.fi„[l], 



+ Kc,fin, (111) 




In dz ^ 



s Jo 




+ ^s?.fi„[l], 



(l-^)Pi 
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The computation of the colhnear-finite parts Kicfinfl] Kc^fin '^i^ be described in more detail in sec- 
tion 16.2.11 and section 16.2.21 

The Vg" configuration (fourth diagram in Fig. I18p shows two additional g-dependent propagators; we 
consider the generalization of Eg. ljlOSp . 



(Pq 



n 



(113) 



and later we will specialize to the case, setting N ~ 2, ki = —q2 —pi, k2 = —qi +P2 and mi = TO2 = Ma- 
We start introducing a Feynman parametrization, 



n 



1 



\ {q + kiY + mf 



T{N) dS^^i{{x}) 



1 



[{q + KY + /l2 



'^{xi^i - Xi) [ml + k1) ~ K"^ , K ^'^{xi-i - Xi)ki, = 1, a;„ = 0. (114) 



Next, we combine the resulting three propagators of Ea. (|113p with variables z and y and perform the q 
integration. Since the configuration in Fig.[TH]is UV finite, we can set n = 4 and obtain 



J^^V{N) dz / dy dS„-ii{x}) 



N-l 



V = [A^yiK^pf] y+m^{l-y), A ^ (K-p z)^+fi'^ . (115) 



Subtracting and adding at the integrand level l/lg", with Vq = Ay + rri^, we can extract the collinear 
logarithm through the y integration. 



J.= -T{N) jdS.^i{{x})j\z I J- (^ln!^-ln^ + |] ^ - j\ 



\ [A-y{K-pf]- 



y 



+ 0(m"), (116) 



where the distribution has been defined in Eq.®. The coefficient of the collinear logarithm can be 
further simplified if we move back to momentum space integrating over all Feynman variables {a;}, 



r 1 " 

r(iV) jdS,,r{{x}) — =Y[ 



1 



^ [ki -pzY + 



(117) 



Formally, wc have reconstructed the product of N propagators on the right-hand side of Eq. (|113p . replacing 
the loop momentum q by —pz. Ea. (|116[ ) has been derived under rather-general assumptions, and it can be 
used for extracting the collinear behavior of any two-loop scalar (respect to q) UV-finitc diagram with one 
external masslcss particle coupled to one light particle; its graphical representation reads 



(118) 



q + p 





= In / dz ( ) + coU. fin. 



(l-z)p 



where the bubble denotes a generic one-loop diagram involving a tensor structure . . g^"* , qa is the loop 
momentum, the wavy line represents the external massless particle {p^ — 0) and the dotted one is the light 
particle with mass m. 
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The result of Eq. (|118p can be immediately applied to the configuration of Fig. [THl note that we do 
not confine ourselves to the scalar configuration, but we consider the full qi -dependent structure appearing 
in the coefficient of the diagram. We obtain 



-p 



s Jo 



(l-z)Pi 




sc.nn 



(119) 



The procedure for dealing with the coUinear-finite part V^"f^^ will be described in section [B. 2. 31 

Finally, let us discuss the extraction of the coUinear logarithms for generic tensor and/or UV-divergent 
two-loop integrals. We introduce 



7" 



iT^^ J {q^ + m'^)[{q + pY + m?] ^J-^ {q + ki)"^ + 



n 



(120) 



and use a Feynman parametrization analogous to the one used for deriving Ea. (|115p . In particular, before 
performing the q integration, we have 



//-"r(iv+2) 



V = q'^ + 2q-IC + M'^, IC = zp + y{K ~p), = {fi'^ + K'^) y + m'^ {1 - y), (121) 

where K and fi^ have been defined in Eq. (|114p . The q integration is performed according to 



/^"'^ n 



T{N + 2-r) 



r{N + 2) 
1 



n 



d 



i—l.r 

r(iV-r + e/2) 



r(iV + 2) 



n 



d 



^-Jv+r-is/2 



(122) 



where V = M"^ - JC'^ ^ [A - y {K - p)^ ] y + {1 - y) is the same of Eq. ljllSp . 

The first derivative in Eq. (jl22p decreases the power of V by one unit and generates an extra factor JC. 
The second derivative acts on both JC and V, producing two terms: one where the power of V does not 
change and the metric tensor appears; another one where the power of V decreases by one unit and an 
additional factor /C is generated. After taking all r derivatives, we find a term containing v^-"-'s/2; 
where all derivatives have acted on V, generating a factor IC^^ . . .IC^''. All other terms in the result will 
contain at least one power of the metric tensor, and for each of them the power of V will be greater than 
—N — e/2, because at least one derivative has not acted on V. 

We can show the following result: if the power of V is greater than —N — e/2 then no coUinear logarithm 
is generated. We cast the y integral of Eq. (|12ip as 



^ 



dy 



{Ay + wP'Y 



dy y" 



1 



1 

^0 



(123) 



where we have added and subtracted in the integrand Vq"", with Vq = Ay + vr? . From Eq. (|12ip . Ea. (|122p 
and the subsequent discussion we can argue that: 1) a < iV -t- e/2 for all terms; 2) /3 > — 1, since 
powers of y in addition to the term of Eq. p2ip can be embedded in the /C factors. Therefore, being 

/3 — a > —1 — e/2 in all cases, the second integral in Eq. (|123p is always finite in the limit m — > 0: indeed 



dy y'" 



1 



^0 



dyy 



[A - y (K - p)^]'^ A°' 



if l3-a>-2. (124) 
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The first term in Ea. p23p . instead, is collinear finite only if /3 — a > —1; in tliis case indeed 

„P 1 



dy 



(125) 



Therefore, for all terms where the power of V is greater or equal to ~N — e/2 + 1, we have f3 — a > — e/2, 
Eq. (|124p and Eq. (|125p simultaneously hold and no collinear logarithm is generated. 

As a result, we can extract the collinear behavior of J^^'"'^'' replacing Eq. (|122p in Eq. (|12ip . and considering 
only the term where all r derivatives have acted on V] in this way we obtain 



(-irr(7V+0(^^j j^dz J^dyJdS,-,{{x})y--' [] y-«-./2 ^ ^^jj (^26) 



The expression for JC given in Ea. (|12ip contains a term proportional to y; it cannot produce a collinear 
logarithm in Eq. (|126p because it increases the power of y and leads to the above-mentioned condition 
/3 — a > — e/2. Therefore, the collinear-divergent piece of J]^^'""'' is simply given by 

JT""- = r (^+ I) f^j dz dy {-zf \{ p""^ y—^l-^ + coll. fin. (127) 

Now we add and subtract Vq " "^^^j integrate in y the added term and absorb in the collinear-finite remainder 
the subtracted one. The y integration gives 



dy ■ 



r(Ar)r(l + e/2) 1 , /m= 
■ in 



e m 

X . 1 - In — 

)! V 2 s 



coU. fin. (128) 



/o (A2/ + m2)"+^/2 r(iV + e/2) A 

Inserting the result in Eq. (|127p and using Ea. (|117p . we finally extract the collinear logarithm for J^^'"'^'' 

-pzY+m 



„ , coU. fin. (129) 



We are now able to compute the coefficient of the collinear logarithm for a two-loop diagram with a generic 
tensor structure in the configuration where only one external massless particle is coupled to a light internal 
one, 




In- 



l--A^(.)--ln- 




dz {-zY p-^i . .p''-- + coU. fin.(130) 



{l-z)p 



It is important to mention the striking correspondence with the special result derived for the simple config- 
uration in Eq. ()106p . 



6.2.1 The collinear-Gnite part ofV^^ 

In order to compute the collinear-finite parts of the diagrams belonging to the V^^ family, we slightly modify 
our approach. We use a complete Feynman parametrization before extracting the collinear logarithms; they 
are then discarded and the collinear-finite remainder is written in a way that allows for direct numerical 
integration. We follow the parametrization procedure described in section 7 of Ref. [47] for the scalar cases; 



38 



the corresponding parametrization for tensor integrals can be found in section 9 of Ref. 
UV poles from these expressions and performing some changes of variables, we get 



] . Extracting the 



'G ^' 



— I dxi dx2 dyi / dy2 (1 — 2:2) 
Jo Jo 

/■I rVi 

/ dxi dx2 dyi / dy2 (l-xs)^^ , , ^ 
Jo Jo 

- / dxi dx2 dyi / dy2 (l-a;2) PG,iPG.j + ^' 
Jo Jo 



h=l 



(131) 



where we have defined auxiliary quantities, Xg — P2 Vi + M'^ and 

X2 (l-a;2) Ag + y2 [x2 Ml + (l-a;2) m\ Ac = (1-yi) [P"^ xi y2 + pi {yi-xiy2) 



'^G,\ ' 
A SC _ 



1 - xi [1 - a;2 (l-t/2)], 



^2 (l-yi), 



1 







In^^ - (lnyi + 2)ln— 
s s 



Ga = l-2;iy2, 
C(2) 



^G,2 



(132) 







2 / dxidx2 / dj/2 3^2 In 



In- 



y2 



W^gIv2=0 



^G,12 - 



^G,22— r, 



1 7 / /" 

- F^i(s) + - + dxi dx2 dyi / dy2 In 
^ o Jo Jo 

nw+^_v,)(l-4)4- 



s 

C(2) ^1 
2 



- 1 2 d, , , ^ 

— in In dxidx2ayi / ay2 — in 



yi ^ 
a;2 



W^G 



W^gU^^o 



In order to get the collinear-finite parts, we have to extract and discard the coUinear logarithms ln(m^/s). 
Concerning the coefficients of the metric tensor, only A^'^22 ^ collinear-divergent behavior; 11 ^^"^ 
^G 12 in a form suited for numerical integration, and there we simply set m = 0. For A^'^22; some of 
the coUinear logarithms are already explicit and we can discard them; the last term of A^'^22 in Eq. p33p . 
instead, requires a further step. Introducing Wg.o = X2 + X2 y2 + y2 m? and Wg.i = ^g + 2/2 M^, we 
have 



f^dx^ 

Jo X2 



In 



Wg 



X2 Wg\x2=0 



^dx2 

X2 
^dX2 
^2 



In 



Wg 

Wg.o 



In- 



W, 



G,0 



W, 



G |3;2=0 



, {l-X2)AG+y2Ml ^. ( Wg,i 
in : lji2 ^ 



Ag + 2/2 M| 



2/2 



1 



. 2 ™^ , , "T-^ 1 2/2 .5 T • 

= — In h In in L12 

2s s Wgi 



A, 



Wg.i 



iln2-^+C(2)+0(m2). (133) 



The integrals in Eq. (ll3ip containing W^^ show a collinear-divergent behavior for to — + at the point X2 = 0; 
therefore, terms where W^^ is multiplied by X2 are not singular and there we can simply set to = 0. If no 
factor X2 is present, instead, we write 









[ dX2 = 


[ dX2 ... -f 


-I 


Jo Wg 


Jo Wg,0 





dX2 I 



/ 1 



1 



^G+2/2M2 



In 



VI^G T^G,0 

^G + 2/2 (Af2+m2) 



WoA 



In- 



In- 



j/2 rn-^ 
Af2 



dx2 

X2 



1 



1 



2 In 



W^G,l 

2/2 



(l-X2)^G+2/2M2 Ag+2/2M2 

-^©(to^). 



+ 0(to2) 
(134) 



Now all coUinear logarithms have been extracted and can be discarded. The finite remainder, however, can 
not be directly integrated numerically, because of the presence of denominators which can vanish inside the 
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integration region. On the other hand, these denominators are always hnear in the variable y2, being both 
Aq and Wg linear in ^2, and all terms generated from W^^ have the following dependence on 



y'i 



y'i 



ay2 + b' ay2 + b 



In 



ay2 + b 



y2 



In 2/2 
ay2 + b' 



(135) 



For them we write 

ln(a y2 + b) 1 



ay2 + b 
lny2 1 



In" 1 



2/2 



In 6 In 1 



2/2 



ay2 + b a 



In 2/2 In ( 1 



2/2 



Li2 (-^2/2) 



1 



ay2 + b 



■ dy. In 1 



■ 2/2 



(136) 



The crucial point in Eq. (jl36p is that in the right-hand side the 1/a is always multiplied by a regulator 
function which goes to zero when the a vanishes. The — iO prescription associated to each mass ensures the 
validity of these expressions and can be used to produce the right imaginary parts of the logarithms. After 
integrating by parts in 2/2, we get for the coUinear-finite parts 



K?,fl„ = 4^c[0, 1,0, 0,1,1], ^.^.finkf] = - Ep5^^s^c[feo,^:4..,4,^,4,^,e;:,„e;,_,j, 



'hk hk 



7;hk 



h,k=l 



where the auxiliary function I^^ depends on six parameters, symbolically denoted by {6}g 



IscMe] = J^dxidx2dyidy2 I j- In 



h 




ymvG\ 


^(■4 


Vg 




Xg ) 








^y.VG\ 




Vg 




Xg J 





2/1 2/2 Pg 

X2 Xg 



2 In 



^Inll 

Pg 



Xc 



21n- 
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2/i2/2??G 

Xg 



yiMl 
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— In 1 + 
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Li2f 


^G 


V Xg 



(137) 



2/i2/2??g \ 
Xg / 



(138) 



We recall that Xg — V2 ?-/i (1^2/i) + 5 in addition, we have introduced the polynomials: 



Pg = xiX2(l-2/i)(P'-p^) - X2M' + r/o = (1 - 2/1) (P' -P2) 



(139) 



Concerning the coefficients of the metric tensor, we have defined A^"]^j^ = A^'^j^j^jm^o; ^^^n — A^ i2lm=0! 
since they are coUinear free; for A^'22 7 start from Eq. (|133p and insert the result of Eq. (|133p . 



.fin 
^G,22 



dxidyi / dy2 



VI 
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-y2Mi 



_l^^2AG±y2Ml 
2 2/2 s 
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Finally, we show the explicit expressions for the arguments of the vector integrals (second equation in 
Eq.|l32l)), 





= -Xl X2 2/1, 
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bla 


= Xl 2/1, 
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p2 
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(141) 



40 



and for those of the tensor integrals (third equation in Ea. (ll37p ). 



= -'2x1 X2V1 {1 - xi [x2 + Vi ^2(1 - X2)]}., h^aM = ^G^i = ^2j/i (1-2:2) (1-yi), 

^1^11 = ^"12 = xiVi (1 -a;2 -X1X2 + 2 xi X2 ?/i 2/2), 

^0^2 = ^g!i2 = -xiX2yi{l-yi), 522^2 = 0, 

^"22 = 2xiyi{l- xiyiy2), b]?22 = = 2:1 yi (1-yi), 

6^2 ^ 0, = [1 - - Xi X2 [l-yi]? - x\ X2 (l-2/l)^ 

= 4^1 = -2:2(l-2;i){l - 2;i[l - (l-2;2)(l-yi)]}, 4^ii = -2:2 (1-2:2) (l-2/l)^ 

4,12 = (l-2;i yi) [1 - 2;i - a;2 (1-yi)], 42^2 = (i-yi) [1 - 2;i - xi X2 (1-yi)], 

4^2 = -2:2 (1-yi) (1 - xi yi), 42^2 = -2^2 (1-yi)^, 

4^22 = (1 - 2:1 yi)2, 42^2 = 4^22 = (1-yi) (1 - 2:1 yi), 

4'22 = (l-^/l)^ dj^n = rflii = 0, 

r/" - - Ti 7/1 Cl -Ti^ rfl^ - rf21 _J22 _ JL2 _ j21 _ j22 _n 

"g,12 — "g,12 — 21 yi 1,1 2,ij, «G,12 — "g,12 — "g,12 — "g,12 — "g,12 — "g,12 — ^! 

^^0^22 = yi (1_- xi yi ^2), ^1^22 = 2;i^yi (2 - 2:1 yi ^2), 

'*^G,22 — '*^G,22 — ^G.22 — '^G,22 — 2^1 2/1 (1 — 2/l)i ^^0,22 — '^G,22 — 

„11 _ pll _ ^1 ^ ^2 12 _ „21 _ „22 _ pl2 _ =21 _ ^22 _ r, 

'^G.ll — ^G,ll — , f^G.ll — '^G.ll — '^G.ll — f^G.ll " '^G.ll " '^G.ll " '-'j 

6^12 = (1-2^1) (l-2;i yi), 4^12 = l-2;i, 

6g,12 = 6^,12 — (1— 2^1) (1— yi)j 6g,12 — ^G,12 — ^' 

4^22 = (1-2^1 yl)^ e"22 = i> 

4^22 = 4^22 (1-2/1) (1-2^1 2/1), 4 22 = 4^22 = l-2/l, 
^G,22 = ^G,22 — (1— 2/1) ■ 
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6.2.2 The coUinear-Bnite part of V^^ 

In this section we consider the V^^^ configuration, using the Feynman parametrization of Ref. [48 . For the 
needed UV-finite scalar configuration we set e = in Eq. (159) of Ref. 48J, perform the y^ integration and 
make the appropriate change of variables in order to simplify the integral. In the coUinear configuration we 
obtain 



/•I ryi 2. 

Kc =- dxi dx2 dyi / dy2 xi (1 - X2) y2 V — — tt— , 
Jo Jo ^^^1^ W^,,U^^^ 



where we have defined 



WK.a = 2:2 (1-2:2) A^c,a +y2 [m2 (1 - ^2) + M2 X2] , 

AK,a = (1-2/1) b? (yi-a;iy2)+pi2:i?;2] + (1-2/2) M2, 

WK,b = WK,a{pi ^ P2), 



UK,a = ^K,a +y2 M2, 
UK,b = UkApi ^ P2)- 



(143) 



(144) 



The integrand behaves as l/x2 for m — 0, revealing the presence of a coUinear divergence. In order to 
extract the coUinear logarithm, we introduce the polynomial WK,aa — X2 A^^a + 2/2 (™^ + Af2 X2) (similarly 
for Wji^b)- Then, we perform the integration over 12, 



' I-X2 

"^2 TJ7 — 







/ ^2:2 




'0 



1 



1 



1 



^x,a+2/2M2 



m2 , M2 



2:2 



-In +ln— 2- + 21n 

s s 



m=0 ^K.a 

A^,a+y2Ml 



m=0' 



2/2 A/2 



A 
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All the polynomials involved in the computation are linear in 1/2; therefore, an additional analytical integra- 
tion is possible. First, we perform two splittings. 
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(146) 
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Secondly, we carry on the y2 integration by means of the relations introduced in Eq. p36p and of the identity 
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Note that for a the content of the squared bracket vanishes; the complete result reads 



sc,nn 



dxidyi 



M2 



L^{yip^,yiMl,yiMl-0) 



2M2 



{pI-pDvi . 

where we have introduced 
L,,{A,B,C;a)^U2(^. 



2M2 



L^{yir],,,0,yiMl,l) 



AP-AP 



Lio 



A + B\ 1 



■ In^ 1 



A 



In^ 



C 



In^ In 1 



A 



(147) 



(148) 



(149) 



2 V Xk, 

Here Xk ~ Vi (l~yi)P2 + -^^^ the polynomials appearing in the arguments of are given by 

rj^ ^L0^ + Ml - AP, = - AP, lu^ = {l-y^){j>l-pl), p^ = {l-y^){AP-Ml). (150) 



6.2.3 The coUinear-Gnite part ofV^^ 

As previously done for and V^^, we specify here the parametrization of the diagram, following the 
procedure outlined in section 10.4 of Ref. [47j . Setting the UV regulator e to zero and performing some 
trivial change of variables, we can cast the coUinear configuration of V^^ as 



where we have defined the following quantities: 



(151) 



WH=y{l-y)AH+yB„ + {l-y)x{l-x)m'^, A^ = x {l-x){l- Z2- z:i){l- zi- Z2) P"^ , 
B„ = xxh + {1-x)AP, Xh=P^Z2{1~Z2) + AII (152) 

= All - A'P ^ {1 - Z„), Z„=zi+x{l-y){l~zi~Z2). 

Since the collinear singularity for m —^ shows up at the point y = 0, we extract the coUinear logarithm 
according to 



1 

dyy 



1 



+ 



1-y 1 



1 



ln^^-ln!^-l 
x{l — x)s s 



f'dyfl-y 1_ 



1 y^+^l-y) /-l 2/"(l-y) , ,0, 2^ 

dy = I dy — — ^0{Tn), n>0, 



+ 0(to'). 
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where we have introduced Wha ~ y y^Hfi + a;(l— a;) to^, = {l—y)AH + Bh and Wnfi — A„+Bh. We notice 
that all polynomials involved in the computation are combinations of the quantities Ah and B^ , both linear 
in zi; therefore, for the terms involving W^q or W^^, we perform an integration by parts in zi using the 
following relations: 



\ „ ln(azi + b) 
dzi 



(azi + 6)2 
^ {azi + bf 



dzi 



{azi + by 



+ n z 



ln(azi + b) 



= — — / dzi 



^ a{azi + b) 
1 



ln(a + b) In b 

—7 — TT- + OruQ —T, 

a(a + b) ao 



azi + b b{a + 6) 



(154) 
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Let now be V and Q two generic polynomials in the Feynman variables. The result after this step can be 
formally written as a sum of terms of the form l/QV or InV/QV. These polynomials are directly related 
to Wh and Wh.o and therefore are linear in Z3, as one can easily prove by direct computation. Therefore, 
we can use the following relations, with A = ad — be: 



1 



1 



{az3 + b){cz3 + d) A 
Injczs + d) _ J_ ^ 
{az3 + b){cz3+d) ' A 



ln(az3 + 6) — ln(cz3 + d) 



, A az3 + b ^. ( A 
In — In \ — ; — L12 



-In- 



a(cz3 + (i) / o(cz3 + (i) 



a cz3+d 

After performing the integration over Z3, the remaining integrands can be cast as 

1 



In 1- 



a{cz3 + d) 



— ln"fn 




1 , In — 






s 



-Li, f ^) 



LOh s 



1 , In a; , ln(l — x) , In — 
s 



.(155) 



(156) 



with n = 1,2, where we have defined loh = x Xh+{}—x) M^. The first two cases in Eq. (|156[) can be numerically 
integrated; the last, instead, requires further manipulations. In particular, for the terms with 1/y, coming 
from the last integral of the first relation contained in Eq. ()153p , we have always the combination 



^r^^ln^:ii^^^-^Li. 

^hJo V V'-LUh UJh 



V 



1 



V'-UJhJ UJh 



Li2 



-In^^ln 1- 



V'-ujh V'-ujh 



V'-UJh 



-C(2) 



(157) 



In the last expression, the terms with C(2) cancel out and the others are suited for numerical integration. 
For the other terms of the type uj^^ \ti{V js), we explicitly write the iO Feynman prescription and get 



I V 1 , {V+ujH)-ujH-iQ 1 
— in — = — In — — 



In 1 



UJh 



r + UJn-iO 



In 



V+ujh 



iO 



(158) 



where the first logarithm regulates the zero of the denominator ujh, while the second is always of the following 
form: 



In 



^i^-iOj =liifx{l-x) — -iO] =lna; + ln(l-a;) + ln 



V 



iO 



(159) 



Here V does not depend on x. As a consequence, for the terms containing this logarithm, we can integrate 
by parts in x, since uJh is linear in x, using the three relations of Eq. (|136p (replacing 2/2 with x). 
Collecting all pieces, we derive the final result for V^^ fin' 



3 . 

V;?,fi„ = E J dC,{x,y,{z}) 



+64 



In 1 



XU), 



2Li2(^) + (ln 

UJ,: 



x(1-x)uJh UJn 



-21n^)lnfl-^ 

UJn 



AP 



In 2 

X s 



In 1 



XUJh 

Jp 



+ {Qs + ee^nx) In 1 



Xh 



(160) 



where, in addition to uJh introduced after Ea. (|156p . we have defined auxiliary quantities in term of Xh of 
Eq.ldSa), 



i^H=XH-M'^, UJn= UJh + P'^x{l-X)an, UJri^y = UJh + P^x{l-x){l-y)an, 



CTl = (1-^2) , CT2 = -Z2 (I-Z2), 



p2 



<^3 = (1~Z2)(1-^1-Z2), 



9„ = 2d„,„ In ( — (I-Z1-Z2) ~iO] + i tt e„,„ sign(P ). 



(161) 
(162) 
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The coefficieiits appearing in Eg. pOOp and Eq. (|162p read as 





= M2-Af2 + £-(4z|-6z2 + l), 


bH.2 


bH,3 


= P^l-X)[l - 2Z2 - 2{1-X){l-Zl-Z2)], 


ChA 


Ch,2 


= p2^[l-f (2-y)(4z|-6z2 + 3)], 


Ch,3 


dn,! 


= PHl-2z,), 

= P^izi~Z2), 


dna 


dH,3 


dnA 


dn.b 


= -P2[l _ 2Z2 - 2X (I-Z1-Z2)], 


dHfi 


Bh,! 


= 4(M2-A42) + i^(4z2_6z2 + l), 

= ^(12z2-18z2 + 5), 

^ P^[x ~ 2{2-x){l-Z2){l~2z2)], 


eH,2 


eff,3 


eff,4 


6^5 


eH.6 



= Ml - M2 - i^[(4z| - 6z2 + 3)x - 2], 

= -2P^x (1-Z2)[1 - 2: (1-Z2)(2~y)], 

= -2P'^x [I - 2z2 - 2 (2-2x+.Ty)(l-zi-Z2)], 

= 0, 

= _ 2z2 - 4(1-.t)(1-zi-Z2)], (163) 
= _ 2z2 - 4a; (I-Z1-Z2)], 

= 2(Af2_M2)-i^(4z|-6z2 + l), 

= -P2[l_2;(422_6z2+3)], 

= p2[^ _ (l-a;)(4z| - 6z2 + 3)]. 



6.3 Vertices with two photons coupled to light fermions 

In this section we extract the coefficients of the collinear logarithms for the Mis illustrated in Fig. [TOl 
where both photons are attached to light-fermion lines. In general the two light masses m and m' are 
different and both vanishing small: therefore, we need to extract both collinear logarithms, L — hi{m? /M"^) 
and L' = ln(m'VM2). 

6. 3. 1 The Master Integral V^^ 

We apply to the parametrization introduced for V" in section 10.4 of Ref. [47]. We first combine the 
qi and q2 propagators with zi and Z2 Feynman parameters; the resulting expression contains the product of 
three propagators, the first in qi, the second in qi — q2 and the last in q2, 

ya _ 1 fn. , P'M^~4{M' + P'+q,-p2)q2-pi 
V^c-J^^jd qid q2j^ dz, dz2 [12]2 [3] [45]2 ' 

[12] = g2 + 2ziqi-pi+m'2, [3] = (gi - (72)' + M^ [45] = ^ 2 gj- (pi + Z2 P2) + ^2 P' + m^. (164) 

The terms [12] and [3] are then combined with a variable x and the qi integration is performed generating 
a new q2 propagator; the resulting two (72-dependent propagators are then combined with a new variable y, 
so that the integration in q2 can be carried on. After extracting the UV simple pole in e, the result simply 
reads 



^cfc ^\fIMI') ~\ - IdC, {x,y,{z}) (1-y) 



ln-^ + (i_2;) ^ 
M2 ^ ' Ug 



Ca = M^[l + 2 Z2 (1-y)] + P' Z2 {\-y) {\ + z^\\-x (1-y)]}, 

Ug = zxz2xv(\~x) (l-y)p2 + xyM2 + a;(l-a;) (1-y) + y (1-a;) m'^ . (165) 

Only the term containing U^^ is singular for a vanishing fermion mass; indeed, for to — > (to' — )■ 0), a term 
y~^{x~^) can be factorized out of 1]^^-, leading to a divergent integral. Therefore, if U^^ is multiplied by 
the product xy, we can simply set the fermion masses to zero, while if U^"^ is multiplied by x or y we have 
to perform a subtraction. 



Ug \Ug UG,y J ^G,a \Ug Ug,x J , 

where the subscripts x and y indicate that inside Ug we have set x^ = and = 0. The first bracket in both 
equations is collinear free while the last term of the first (second) formula generates the collinear logarithm 
L (P'), after an explicit integration in y (x). For terms where the coefficient of [/^^ contains neither x nor 
y a double subtraction is needed, 

1/1 1 1 1\/1 1\/1 1\ 1 

'^G V'-^G '-^G.x '^G.y '~'G,xyJ \L^G,x ^G,xyJ \'~JG,y '~'G.xyJ '~'G,xy 
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with — — m the polynomial with the subscript xy. The first bracket of Eq. (jl67p contains no 
singularity; in the second (third) bracket, we integrate in x (y) generating the collinear logarithm L' (L). 
Finally, in the last term, both integrations in x and y are performed and the product L L' appears. Having 
extracted the collinear logarithms, we can now set the fermion masses to zero, obtaining 

Vl = ^Fl,{M^) + ^+ j dC,{x,y,{z}) {laLL' + IlL + IlL' + (168) 



where we have introduced 



1 i„ ^G,0 



Co In 



M2 



■jxy 



xy 



++ 



izluM^G.uiM^) - 1] 



^G,0 



^G,0 

M2 



u = x,y, 



In" 



C(2) 



SG.a; 



1 ln(l - x) 



(2 In " 1) ~ CgI ln(l ~x) ~ ac^J - {l-y)\n^. 



(169) 



u—x^y 



Here the '+' and ' + +' distributions have been defined in Eq.® and Eq.®; the coefficients and the 
functions are given by 

aa^P^Z2{{2~y)[l + zi(2-a:)] + z^{l~x){l-yf] + A'P[l + 2z2(2-y)] , S^^^ z^Z2{l-x){\-y) + 
ctG,y = ctG{x = l), q;g,x = Q;g(2/ = 1), aG,i = "G(a;=l = y), 

?G,, =?G(a: = 0), ^G,x =fG(y = 0), ^Gfi = ^G{x = Q = y). (170) 

Since the ^g functions are linear in the four integration variables, all integrations can be analytically per- 
formed. Introducing u) — —P'^/M'^ and = ln(l — oj), the result is 



LL' [2w+ (1-w) (2-w)U + L [2w(l+cj) + 2 (l-w)L - 



+u;(l-w)Li2 (w)] +L' [6u;+ (1-cj) (4-cj)L - (1-^^) (2-w) + (w2-2cj + 2) Lis (w) " 



19 

f +4C(2) 



-(5+c.) + 2(2-c.)C(2) 



(1-Cc;) L - (1-cc) (6-Cc>) /2 + - (1-c^) (2-c^) ll 



" (tj2-4cj + 2) Lis (t^) + 2 (1-t^) (2-t^) LLis (w) - 2cjLi3 (cj) - 2 (3tJ-2) 5i,2 (t^) }>, 
where Li„ (w) and 5*1. 2 (ti^) are Nielsen poly- logarithms. 



(171) 



6.3.2 The Master Integral Vf^ 

We parametrize Vf^ following section 10.4 of Ref. [17]. After combining the qi, qi — 52 and 52 propagators 
with Feynman parameters zi , 22 and Z3 , we obtain 



2e 



d"gid"g2 / dCs {{z}) 



p2Af2 + 2p2gi -pi -4(gi -^1)2 



[12]2 [34]2 [56]2 



[12] = ql- 2 zi qi ■ p2 



[34] = (gi - 92)' + 2 (qi - (72) • {Pi -Z2P) + M2 



g| - 2 (I-Z3) 92 -Pi + n^'^- 



(172) 



Afterwards, [12] and [34] are combined with an additional variable x] we integrate in qi and obtain a new 
92 propagator which is combined with [56] using a variable y. Finally, the 52 integration is performed. The 
MI V^" is UV finite; after setting e = 0, we obtain 



= -/ dC^{x,y,{z}) xy{l-x){l~y)^, 
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c„ - P2 ^^2 ^ p4 _ a; (Z1-Z2)] [1 - zi + a; (1-y) (zi-za)], 

i/^, = (l-x) (zi-Z2)(l-Z2-Z3) P^ + xyx + y (l-x) + x {l-x) {l~y) m'^ , (173) 

where we have introduced % = Z2 (1 — 2:2) + A'P . 

The colhnear-singular behavior corresponding to to — > (m' — > 0) shows up at x = (y = 0). In fact, 
in this hmit, the global xy factor combines with the l/x^ {l/y"^) term coming from U~'^ and produces a 
global x^^ (y^^) divergent integral. Subtraction terms are evidently not needed when the coefficient of U^'^ 
contains the product x^y^; for terms where J7~^ is multiplied by either x^y, xy^ or xy, we apply the following 
subtraction formulae: 



x^ y (l-x) (1-y) ^ ^2 y ^-2 _ ^-2^ ^^2y ^^2^ 

xy" (l-x) (1-y) = .^y2 _ ^-2^ ^^y2 

xy (l-x) (1-y) U-^ = xy [(l-x) (1-y) C/-^ - (l-y) U'^l - {l-x) U'^^ + C/^^^ 

r-2 



-2 
-2 



+ X y [(1 - y) C/^j - Ujy] + X y [(1 - x) C/^j - + x y 



(174) 
(175) 

(176) 



The subscripts x, y and xy for Uh indicate that we have set x^ = 0, y^ = and x^ 



respectively. 



The first squared brackets in all three Eqs. (|174m76p are coUinear free. In the second (third) squared 
bracket of Eq. (|176|) and in the last term of Eq. (|175|) (Eq. (|174p ). the integration in x (y) generates the 
coUinear logarithms L {L'). In the last term of Eq. p76p we perform both integrations in x and y obtaining 
the product LL'. At this point, the coUinear-divergent behavior is explicit and we can set m = to' = 
inside all polynomials, getting 



Vf. = J dC, (x, y, {z}) {I^LL' +IIL +IIL')+ F-.fi„> 



Here V^" reads 



f-2 



x,y. 



(177) 
(178) 



E 



Co(l-^||-2 1n%| + C(2) 



1 



M2 



(l-x)(l-y)^-^ 



xy 
2ln(l-x) 



++ 



, (179) 



where the an and Ph coefhcients and the functions are given by 

a„ ==P4 (l-x) (1-21-22) {(1-21-22) [1-2: (2-y)]-zi + Z2}, /3« = [P^ 21 (l-zi) + M^] , 
e„ = P2 [(1-^1) (I-Z2) (I-23) + 21 22 23] (l-.x) (1-y) + P2 22 (1-22) (x + y - xy) + M^, 
aH,y = anix^l) =^ 0, Q;„,a; = Q;„(y=l), ^^^^j, = ^„(x = 0), Caa: =CH(y = 0), ^„_o = CH(a; = = y). (180) 

We consider now the coefhcients Xh, 1% and X^, written in Eq. ()178|) . of the coUinear logarithms of Eq. (jl77|) . 
and express them in terms of one- loop functions. After changing y ^ x in Xj^, we observe that the coefficients 
of L and L' are integrals in x, 21, 22 and 23. Then, we use the following trick to get rid of ln(^„.o/-^^)i 



P2 22 (1-22) +M2 



\ (i-^)eHi 

ax — 



(181) 



Since Ph and aH,x do not depend on 23 and (^h,x is linear in 23, the 23 integration can be easily performed, 
leading to 

Vi', = fdx dzi dz2 ( Jf LL' + JIL + Jl' L' ) + jf = -/3„77->-io, 

Jo 
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1 P^Z2{1-Z2)+AP 



V„.oPh.o In 



M2 



, Ji" = + (182) 



where we have defined 

Vh.x = P^Zi Z2 + P2 (1_22) a; + rj„.o = P^^l Z2 + ]VP, 

Ph,x = (I-Z2) (l-x) + P2 ^2 (I-Z2) a; + M2, p„,o = (l-^a) + 

At this point, we perform some partial fractioning, 

Ph (1-2:) Vh!x Ph% = {{zi -XZ2) r]-]^ + [1 - zi - X (I-Z2)] Ph^x}^ 
PhV^I Ph,o = [^1 V^}) + (l-^i) Ph,o] > 



ttaa: 77£f,J; Ph.o; 



= [(I-Z1-Z2) (1 - - zi + Z2] (?/^^|^ - 



(183) 



(184) 



Next, we change zi ^ 1 — zi, Z2 ^ 1 — Z2 in all terms containing p^^^ or Ph,o, which will then become 77^,0 
or 77ff,Oi and obtain 



X 



4P2 Z2ry;i + 2P2 zi ,7-1 In 



.2. „-i ,„P'^2(l-^2) + M2 



M2 

--2p2zi,7-i. 



J-/- J^ + 2p2[(l-zi-Z2)(l-a;)-zi + Z2]r;-_i, ^i-^ ' - -2 P^ zi 77-^. (185) 

Finally, we can get rid of the logarithm in Eq. (|185p by using a trick similar to the one introduced in Eq. (|18ip . 



_1 P2 Z2(l-Z2)+Af^ 

^H,o In Jf2 



1 „-i 

, 'IH.X 

ax — — 



, -1 , ^H,0 



(186) 



For Z2 rjjj \ in the first line of Eq. (|185p we use instead 



2 P2 Z2 = P2 (1-Zi) 77-1 + P2 Zi ^ 



+ X 



-d^. In . 



(187) 



In this way the terms containing the '+' distribution cancel in Eq. (jl85p . In all terms which contain ri^\.^ we 
first change Z2 —* Z2/X and then zi — )■ z, Z2 — > X2 and x ~* 1 — xi + X2- In the rest we are able to integrate 
in all variables x, zi and Z2, obtaining 



1^." = 2 1 



1+w 



L, LL' 



1 + L {L - 1) + Li2 {u}) 

UJ 



{L + L') 



rl rxi 

-2P'^ dzdxi dx2{{l-z)L+ [l-a;2 z-a;i (1-z)] L'} + ^^0,(111: 
Jo Jo 



(188) 



where we have introduced oj = —P'^/A'P, = ln(l — lj) and = P^ a;2 [1^X2 z—xi (1— z)] + Af^ (1— a;i+a;2); 
the latter can be easily recognized as the polynomial of a one-loop triangle with z-dependent momenta. 
Therefore, we can express all coefficients of the collinear logarithms in terms of integrals of one-loop functions 
as for the single-collinear case treated in section |6^ 



Vf = 2[l-^L\ LL' + 2 



1 + L {L - 1) + Li2 (uj) 

UJ 



{L + L') 



P-2 + 2Pl 



2 fdz [{l~z)P^L+{P^ + 2q-p2)L'] 
Jo 




(l-2)Pl 



* ac.nn 



(189) 
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In Ea. (jl89p the triangle of loop momentum q is to be intended in the graphical representation given in 
appendix [X] 

For the coUinear-finite part V^^ we do not need an analytical expression, since it will be enough to 
put it in a form suited for numerical evaluation. The situation is similar to the one for VJ^f^^^ computed in 
section [6.2.31 all terms of Ea. (|179p have the following structure, with respect to zi and z^: 



In^ ^ 

f{n,m) = z^ ^2~' 



— a zi z^ + b zi + c z^ + d, 



n,m — 0, 1, 2. 



(190) 



Therefore, we can integrate by parts in zi after inserting 

2 



e az3 + b 

obtaining the result 



hi 
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^ -d.,l, (191) 
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where we have used the following short-hand notations for the ^ functions, 



(192) 



C = Cff, S,x=iH,x, iy^^H.y, Co^CffiO, 6 = -P^ (^^2 + ^3 " 1), 

CA = e(^i-o), = ex(2i = 0), = ey(2i = 0), eoA = eo(2i -0), (193) 

Cb=C(zi = 1), C:eb = Ci;(2l = 1), Cys = ^2/(^1 = 1), Cob = Co(^;i = 1), 



and for the a and (3 coefficients. 



a, = p4^i{l-2(l-^i-Z2)[l-x(2-2/)]-2(l-zi)}, ai = ^2 [l - a: Z2 (2-y)] , 
Oixz = [1 - 2 [l-x) {I-Z1-Z2) - 2 {l-zi)] , = P4 ^2 (1 - XZ2), 
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(194) 



a.o = ^'ni-^2)[l-a:(l-Z2)], 



/3, = p4^_2zi), /3o=p2Af2 



The result contains terms having the following structures in Z3: l/QV and 1 / QV In" V, with n = 1,2, where 
Q and are linear in Z3. The integration in Z3 can be performed using Eq. (|155p for the first two cases and 
the following relation for the last one, with A = ad — be: 



\n^{cz3+d) 
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{az3 + b){cz3+d) A^^n^"^ a 



-21n(cz3+d) Li2 



ln(az3 + fe) — ln(cz3 + (i) 
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2Li3 
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a{cz3+d) 
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ln^(cz3 + (i) 



In 1 
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a{cz3+d) 



(195) 



We also integrate analytically the '+' distributions and use the properties of the dilogarithms to cast the 
result in one of the following forms: 

:^ In'^fl-^^") , n = 1,2; -Lisf-"); -uJ-); -Sx2(-); - In ^ In'^^, fc = 0, 1, 2, (196) 
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where is a polynomial or a logarithm of polynomials. The first four forms are ready for numerical 
integration, since the zeros of x are compensated by a regulator function. For the last case, we observe that 
the polynomial V is always of a special form: making explicit the i Feynman prescription we have 



■In- 



M2 



In 1 



X 



V + x-iO, 



In 



V+X 
Af2 



V + x = ±aP^ (I-Z1-Z2) 



0,1, (197) 



where a is the product of non-negative factors (1 — x), (l — y), z^ and (1 — 2:2) which can be extracted from 
the logarithm. Our procedure has generated a regulator function for the zeros of x a-iid obtain a simple 
remainder which can be integrated in one variable (in our case z\ ) . 

Let us consider now terms without a regulator function: thanks to the special form of V^j^, the integration 
generates a factor x in the numerator, which cancels the x in denominator. Note that this cancellation 
was anyway expected: if V^^ is a master integral, then it must have the right threshold properties of the 
original diagram. From the beginning (see section [5|), we do not expect threshold singularities in and a 
non-regulated 1/x factor, having a l/-\/l — AIvP / s threshold behavior, is not allowed. The final expression 
for V^" reads then 
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(198) 



where the polynomial in the denominator is x 
are given by 



Z2 {I — Z2) P'^ + A/2. The ^ functions and the /3 coefficients 
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(199) 



/3a,i=4[l-2(2-2a; + a;t/)(l- 

/3a. 2= -4z2(l-2a;z2 + a;t/2;2), 
/3a,3=2-a:(l-2z2+4z2)(l^ 



zi~Z2)], /3fc,i=4[l-2(l-zi)-2(l-a^)(l-^i-^2)], Pc,i 
/3h,2= -4^2 (1-2^2:2), /3c,2 
-y), (3b,3='2-x + 2xz2{l-2z2), /3c,3 



l-2zi, 
-A/2 /p2 (200) 

Aiyp^. 



6.4 The special case of the family 

The computation of the amplitude for // ^ 77 has been performed keeping the light-fermion masses 
during the whole generation (g) simplification procedure; they have been set to zero just after the extraction 
of the coUinear logarithms. Therefore, at intermediate steps, we have diagrams multiplied by a small mass 
which, in general, will disappear when the fermion mass goes to zero. However, some of these diagrams belong 
to the V" family, and develop a l/m2 behavior. They are shown in Fig. and they would in general lead to 
a residual contribution when their coefficient is proportional to m^. It is worth noting that, concerning the 
H 77 amplitude, only the scalar V^' function remains after the reduction® symmetrization procedure; 
in fact, all tensor structures have been reduced to simpler topologies. 

We parametrize the T^" {i = a,b) integrals shown in Fig. [20l following section 8 of Ref. [47j and obtain 

= - (^^^ r (2 + e) J\xdy J^zi J\z2[x{l - x)]"! (1 - y)^-\y - z,)x^{x, y, 21,22)-'"% 
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Figure 20: The V^' family, characterized by a 1/m^ behavior in the m ^ limit. Dot-lines have a small mass m; 
wavy lines are massless. 
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with TO^ ~ [-^^1(1 — x) + M|x]/[a;(l — x)]. In order to extract the aforementioned 1/m? behavior, we apply 
a Bernstein-Sato-Tkachov algorithm [84 on Xi, 



X^ix,y,zi,Z2y^^^ = ^ <! 1 



1 



1 + e 
1 

1 + e 



1 



1 



{y - l)dy + -zidz, + -Z2d, 



xAx,y,zi,z2) ^ % 
X,(x,2;,Zl,Z2)-l-^ (202) 



After integration by parts the poles are extracted. Analyzing the behavior at small m, we find that, in 
addition to the 1/m? behavior, also Inm^ terms are present. The structure of the remaining integrals is 
simple enough to allow for a direct integration in all variables. The result yields 



X — -^'^ 



A/2 m ' 



1 

P2 



2F^,(.)-ln(^-^j +p(-P^Mf,M|) + l 
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771 X www:: 
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P2 



[l + m- p{-pI Ml Mi)] -2-in + 



(2 + ITT) p(-pi Mf, M|) + 0(e) + 0(m2), 



-Pi, 



(204) 



and with p(x, y, z) — [y hi(jj/x) — z ln(2/x)]/(2: — y\ Therefore, before taking the massless limit in the terms 
of the amplitude generated by the V^' family, we introduce the results of Eq. (|203p and Eq. (|204p . 

It turns out that all terms coming from these kinds of V"' diagrams, where the Higgs boson and the two 
photons are coupled to light fcrmions, are either zero or cancel analytically one another. As expected, they 
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lead to a total vanishing contribution; we would have obtained the same result treating them as massless 
diagrams from the very beginning. Note that diagrams with leptons, up and down quarks of the first and 
second generations and bottom quarks cancel separately. 

Summarizing, in this section we have discussed an algorithm which allows to express the coefficients of 
all coUinear logarithms in terms of one-loop integrals. In order to verify their cancellation, we then reduce 
all one-loop diagrams to scalar integrals with propagators raised to a canonical power one; this last step 
requires the standard reduction of Ref. [FT] and the IBP identities of Ref. 41j. Each coUinear-divergent 
diagram contributes to the total result with its coUinear-free part which is manipulated as described above 
and then computed numerically. 



7 Evaluation of massive diagrams 

For the evaluation of massive two-loop diagrams, we follow the methods of Refs. |45l46j for self-energies 
and of Refs. [47148149] for vertices. For the three vertex families V^, and V" , we modify our approach, 
since two external massless particles appear and simpler results can be derived. Since the massive diagrams 
do not lead to coUinear logarithms, the techniques described in this section can also be applied in the case 
of external light particles (setting their squared momenta to zero). 



family. Let us consider the scalar, vector and tensor integrals of the family, in the configuration 
where coUinear divergences are absent and one of the two external squared momenta and p| vanishes. 
Setting without loss of generality p| = 0, we follow the parametrizations described in section 7 of Ref. [17] 
and in section 9 of Ref. [IS] , in close analogy with the procedure already described in section 16.2.11 The 
extraction of the UV poles is trivial; after an appropriate change of variables, we can write the functions 
in parametric space as 

/•I rvi 
V^IpI^o^ - / dxidx2dyi / dy2{l-X2)Yj'^, 
Jo Jo 

VWi]pl=o = / dxidx2dyi / dy2 {1~X2) p'^^,Y-\ = Va^^iP^, 

Jo Jo 

V''{l'^qj]pl=o ^~ [dxi dx2 dyi / d2/2 {l-X2)p'^,pl., + ^a.j (205) 
Jo Jo 



where Yq reads 

YG^aayi + ba, ac = (I-.T2) ctc (Tg = X2[xiy2{pi - P'^) + m\- ml], 

bo = xi X2 y2il~X2)[P^-pl{l-xi + y2-X2y2)] + a;2(l-a;2) (7^5-77142/2) + 2;2(l-a;2)xG.x + "13X27/2, 

Xg,x = pIxi{1-xi) + m\xi + 7712(1-2:1), (206) 

the coefficients of the tensor structures involving the external momenta are given by 

0^,1 = 1 - a;i [1 - 0:2 (1-^2)], a|a = a:^2 (1-yi), 0^,2 = 1-2^1^2, a| 2 = l-2/i, (207) 
and the coefficients of the Kronecker delta are 



1 TV? 1 f^^ ( Y I — X2 Y 

- / dyi{lnyi + 2) In + - dxidx2dyi / dy2 X2 In — H In ■ ^ 

Uo -5 27o Jo \ s t/2 

1 7 1 Y 

F-i{s) + -r7 + 7; I dxi dx2 dyi / dy2 In ■ ^ 
16 2 Jq Jo 
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lF\is) + lF\is)(^J^'dx,\n^ 
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s 4 
yi 



---cm-- 

32 8 ^ 4 



dxi In 



-- / dxi In ^^£i + t; I d,xi dx2 dyi I dy2 — In ■ 



Ya 



X2 Yc 



(208) 



with = TO4 yi + mi,{\~y\) and s = . 

The terms proportional to 8^'^ in Eq. (l208p can be directly integrated numerically; for those terms in 
Eq. (|205|l containing we observe that the polynomial Ya is linear in y\ and integrate by parts, 



In 1 + -^ -y^Mn 1 



ri / dyi In ( 1 



^G 



Therefore, the expressions of Eq. (|205p can be written in terms of simpler integral representations, 



(209) 
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/I 2 Mi^i2 1 u \ 
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In 1 



aG2/2 
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-aG,zaGjP2P2 



2 2/i(l-2/i) -^j^/"-^ , aGyi\ (1-2/2) 



Og 



■Inl 



ctg 



aG2y2 
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where we have introduced short-hand notations related to Eq. (|2G6[) and Eq. (|2G7p . 



*G — "Gla2^aia2J 



be — ^G|j/2-*yiJ/2J 



^G,2 — '^G,2j 



«l.2=l- (211) 



Looking at Eq. (|210p we see that the zeros of ctg in the denominator are smoothly compensated by logarithms 
that vanish for CTg = {ug is proportional to ctg)- However, in the process under consideration we have 
TO4 = TOs, and we can encounter some numerical instability; in this case, CTg is proportional to X2 2/2 and 6^ 
vanishes for X2 = 2/2 = 0. Since the zero of &g is of lower order with respect to the zero of CTg, the instability 
can be cured by a sector decomposition of the unit square [85| . 



Jo 



Jo 



J/2 



dx2dy2 f{x2,y2) ^ { dx2 dy2+ dy2 dx2 ] f{x2,y2) ^ / da; 2^2/2 a;2/(a;2, a;2y2)+ 2/2/ (2:2^2, 2/2) • (212) 



family. We analyze here the configurations with — and Pi=P2 — 0, focusing on the 



scalar diagram, since it is the only one which survives after the reduction procedure described in section |4?T1 
Note that the method explained here can be applied, with a straightforward generalization, also to tensor 
integrals. We set e = in Eq. (159) of Ref. [48j, perform the ya integration and make some change of 
variables simplifying the integral. 



/•I ryi ^ 

/ dxi dx2 dyi / dy2 Xi X2 
Jo Jo - „ , l^K,i 



X2 (l-2;2) 



1 



a;2(l-a;2)^K,i + 2/2Sjf,i -4jr,j 



(213) 
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AK,a = P^xiy2{yi-xiy2) + {yi-xiy2)[plil-yi) +ml] + xiy2[plil-yi) + ml] + il-yi)ml, 
BK,a = P'^x\xi{l-xi) + X2[m\{l-xi) + m^xi] + {l-X2)m\ ~ X2{l-X2)[m\{l-xi) + m^xi], 
AK,b = AkA'^i ^ m2,m4 ^ m6,pi ^^2), 'Bk,;, = SK-,a("^i ^ m2,m4 ^ m6,pi ^^2)- (214) 

Note that the zeros of BK,a and of S^.b in Ea. (|213p are compensated by a corresponding zero in the squared 
bracket. The general procedure described in Ref. [47] writes the integral in yi and y2 as a one- loop C 
function; however, after setting — or p\ = p^ = 0, we observe that Akm and A^^b are linear in 2/2 and 
2/1, respectively, thus allowing for an integration by parts where we use the last relation of Eq. (jl36p . For the 
case = we obtain 



dxidx2dyi^^ 
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ml IH yi 
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^Pl ^P2 J 



«A-,a = b2--Pi)(l-yi) + - "^4> ^K^a P?2;i(l-2/i) + m4yi + TO5(l-yi), 
A similar result holds for the configuration where Pi ^ P2 = 0, 



X2(l-a;2). (215) 
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(216) 
(217) 



X2(l-a;2) ' 

There are special configurations where the result can be further simplified integrating in one additional 
variable. This step becomes necessary to achieve numerical stability in all cases where special relations 
between masses move some zeros of the denominator to the border of the integration domain, as in Eq. ()216|) 
when 7713 — 0, mi — and m2 = mg; in this case we can factorize a;2 from BK,i and carry on the integration 
in 2:2, 



dX2 

X2 



In 1 



AX2 

1-X2 
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BX2 

I-X2 



-Li2(l-yl) + Li2(l-B) 



(218) 



family. For evaluating diagrams of the V" family without coUinear singularities, we rely on the 
methods developed in Ref. [3] where the requirement is that two external squared momenta vanish, as in the 
process H — > 77; this condition allows for a double integration by parts in parametric space. However, the 
method described in Ref. ^ needs special care when applied in the region above threshold (a case not covered 
in Ref. [3])- Following the same parametrization procedure described in section [6.3.21 we set pi = P2 = 
and obtain for diagrams up to rank two 



where is given by 



dCs {x,y,{z}) xy{l-x) {l-y)Y^ 



/ dCs (x,?;,{z}) xy(l-x) (l-2/)y^2p^ ,, ^^a^P/^' 

h=i 



dCs {x,y, {z}) 



i-y 

Yh 



xy{\-x)Y~^p^^^^vl^.^ + -a\^^8^^ 



(219) 



Yh^ xy{\-x)(\-y)xH,^ +y(l-a;)x^i ^VXh.i + a;(l-a;)(l-y)x„,3, Xh,o = {\-Z2-zz){\-zx-Z2)P'^ 

Xh,i = '71i(1-Zi) +7772^1, = ^2 (1 - ^2) + TO3 (1 - 22) + TO4 Z2, Xjf,3 = "^5 (1 " ^s) + TOg^S , (220) 

and the coefficients of the tensor structures read as 



'-ff,i 



^2;(1-2;)(1-Z2-Z3), 



As = 1-zi-a; (l-y)(l-Z2 -zi), aj,_2 = Z3+y (I-Z2-23), 
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'■h,2 



y (I-Z1-Z2), 



7° 



'■H.12 



xy, 



''H.22 



(221) 



Focusing on the configurations which appear for the — s- 77 decay, we set m2 — mi, — m^, niQ = 
and use the symmetries of the result; vector integrals can be simply written in terms of the scalar one, 



1 



1 



where the subscript s denotes the configuration = p\ = 0, m2 
the expression for shows a simple dependence on zi and Z3, 

= [zi Z3 - (1-Z2)(zi + Z3) + (1-22)^] + fin, 
Ph = y{l-x)ml +xyxH,2 + ^ {^-x){l~y)ml, 
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TOi, m4 



ms, me 



(222) 

rris- Furthermore, 



Ih = xy (l-x)(l-y)P^, 
ff.2 =P^^;2 (1-2:2) +m3. 



(223) 



This structure corresponds to the coUinear case VJ^ discussed in section [^231 integrating by parts in zi and 
Z3 according to Ea. (|154p and Eq. (|155p . we obtain (n, fc 7^ 0) 
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(224) 



where in the last equation, which is related to the tensor form factors proportional to 6'^'^ , we can carry on 
numerical integration. For the other four expressions in Eq. (|224p we have to pay attention to the zeros of 
Ph (the factor cancels with the numerator in Eq. (|219p ). It can be easily proved that the real part of 
the expression contained in the squared bracket vanishes in the limit P^ — > 0; however, an imaginary part 
can in general survive. Below threshold, for s = —P^ < 4m§, Ph does not change sign in the integration 
region; therefore, for the real part of the diagram below threshold (as done in [3|) it is enough to take the 
real part of the integrand and then integrate numerically. For H —^ 77, we are interested also in the behavior 
above threshold, where we have to take into account the imaginary part of the squared bracket. Following 
the approach of section I6.2.3[ we make explicit the i prescription for masses and perform the following 
symbolic decomposition: 
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The first two terms are regular when — > 0; for the last term of Ea. (|225p we notice that Ph is linearly 
dependent on y and in the argument of the logarithm a factor y can be simplified between Qi and Q2- 
Therefore, we can integrate by parts in y getting 
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with ui ~ 1 and U2 — y. Here we have introduced short-hand notations for the [3 factors, 

fJl = (1-Z2)^, (72 = -22(1-2:2), 0-3 = (l-Z2)(l-2l-Z2), 0-4 = (l-Z3-Z2)(l--Zl-22), (227) 

for the O -functions and for the a, (3, 7 and J coefficients 

e"% = l3'^:!; + iif,e{zi+z2-i) + 6';ifje{z^+z2-i)e{z3+z2-i), 1 = 1,2, (228) 



aj;\\ = 2x2(l-y)2(l-z2)2, = -^(l-y)2(4z|-6z2 + 3), 

<\ = -4x2(l-y)2(l-zi-z2), ^ 0^ 

««4i="H,n= -^(l-2^)(l-^2)[l-2a:(l-y)(l-Z2)], = a^'A = -f(l-y)[l + ^(l-?/)(4^|-622 + l)], 

"ff,n= "l.n== a;(l-y){l + 2(1-Z2)[1 - 2x{l-y)]}, = = -a;(l-y)[l - a;(l-y)], 

= l-2x(l-y)(l-Z2)[l-x(l-2/)(l-Z2)], a^n - -[l-a:(l-y) + ^(l-t/)2(4z|-622 + 3)], 

a^22 =_2[l_a:(l-y)][l-2zi-2a:(l-y)(l-zi-22)], <A = 0, 

(229) 



/3i,n = 0, = 2x2(l-y)(l-Z2)(l - 2z2), 

7in = 0, 7^," = 2x2(l-y)(5 - 4zi - 4^2), 

'^in = -^l^i = /3i,\2 = = 0, /?|\^ = Pl^^ = f [l + 2x(l-y)(4z|-622 + l)], 

7^12 7^,21 = 0. 7|\2 = 7ff,2i = -a;(l - 2z2)[l-4a;(l-y)], 

S]ii2 = 5)1^1 - 0, <S|V2 = J^^A = -2a:[l-2x(l-2/)], 

/3i,22 = i /3I22 = - 2a;z2(l-y)(3 - 2^2)], 

7^,22 = 3-4zi, 71^22 = -2x[5-4zi-222-x(l-y)(7-40i-4z2)], 

'5l;,22 = 0. 4,22 = 0- 

(230) 

Note that we have presented here only the coefhcients of V^"[q'^q'[\, the surviving configuration after the 
reduction ® symmetrization procedure. Other tensor structures can be obtained with the help of the same 
procedure. 



8 Numerical results 



In this section we discuss the numerical results for the two-loop corrections to the decay widths of the 
processes H gg and _ff — > 77, given by 



T{H ^ gg) = 



4\/27r3A/„ 
1/2 



AQg 

^phys 



and T{H ^ 77) = 



32V27r3A'f„ 
1/2 
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(231) 



^ith A^^y^ ( V2 Ml ) = A from Eq. m and AH^^ ^^(V^G, M^, ) = A from Eq. m ■ 

One important production mechanism of the Standard Model Higgs boson at the LHC is the gluon-fusion 
X ^ H + X . Its partonic cross section ct, to LO in QCD, can by related to ^p^ys 



channel, pp ^ gg 



<99 ^ H) 



32V27rA//4 



A99 
^phys 



(232) 



The relative correction 5, induced by the higher order corrections, is given by F = ro(l-l-(5)(cr = cro(l-l-5)), 
where Eg (ctq) is the lowest order quantity. For the decay H — > 77 we split the relative correction in 5 = 
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(5ew + ^QCD, to distinguish the contributions arising from electroweak and QCD corrections. For the process 
H gg{gg H) we consider only two-loop electroweak corrections; interplay between electroweak and 
QCD corrections has been discussed in Ref. [55| where an estimate of the remaining theoretical uncertainty 
is also presented. 

As discussed in section [5. II and section [5.21 threshold singularities appear in both reactions if we restrict 
the calculation to the RM scheme (see section [573]) . For the case iJ ^ 77 we consider a phenomenologically 
relevant Higgs-mass range from 100 GeV to 170 GeV, which contains the WW threshold. For the gluon-gluon 
case we consider a Higgs-mass range from 100 GeV to 500 GeV. Here we cross not only the WW threshold, 
but also the ZZ and tt thresholds. The WW and ZZ thresholds have a different behavior as compared to 
the tt threshold. Indeed, the amplitude around the tt threshold contains a potentially dangerous two-point 
function, as explained in section [5l but in contrast to the WW and ZZ cases this singular function (for 
= 4 Mf) is here protected by a multiplicative factor = 1 — AM^ /s. In order to cure these singularities 
we have introduced complex masses [60] as described in section 15.31 We stress the fact that at the single W 
and Z thresholds no special enhancement occurs. 

All light-fermion masses have been set to zero in the collinear-free amplitude and we have defined the 
VF-and Z-boson experimental complex poles by 

s, = ifi, -ijj), /i ■ = - r,^ jj = r, (1 - r2/(2 Mf )) , (233) 

with j — {W, Z}. As input parameters for the numerical evaluation we have used the following values [86|87j : 

M„ = 80.398 GeV, = 91.1876 GeV, Mt = 170.9 GeV, r„,- = 2.093 GeV, 

= 1.16637 X 10-5 GeV~^ a(0) = 1/137.0359911, as (M^) = 0.118, = 2.4952 GeV. 

For r{H gg) and cr{gg H) the behavior of (Jew as a function of is given in Fig. [2TJ For the 
gluon-gluon case, treated in the CM scheme, we observe a smooth behavior in the full range of where all 
cusps present in the MCM scheme have disappeared; we may conclude that results in the CM scheme nicely 
interpolate those of the RM scheme around thresholds. Therefore, all pathological aspects associated to the 
crossing of thresholds have disappeared. It is also worth mentioning that the CM scheme greatly improves 
stability in the numerical evaluation of all master integrals, especially those in the V" family. 

In the shown Higgs-mass range the percentage correction varies between about -)-6% and -4%. 

It is important to consider this result in more details: around the WW threshold we find a maximum for 
the total electroweak percentage correction, and only a light shoulder of the curve around the ZZ threshold. 
Both characteristics are less suppressed if one considers only light fcrmions with real masses, shown in 
Fig. [531 Let us define light-fermion correction as those coming from the first and second generation of light 
(massless) quarks and the one from bottom quarks only; each of these three classes of diagrams constitutes 
a gauge- invariant subset. 

It is worth noting that below the WW threshold the contributions from top quarks are small but they 
become significant above the ZZ threshold and even more important around the tt threshold, where the 
curve exhibits a minimum, which is absent for light fermions, see Fig. 1231 In this region the top-quark 
contribution leads to a sizable effect of about -4%. 

The results for T{H gg) have been used to derive NLO electroweak corrections to Higgs production 
at hadron colliders in Ref. 88]. 

The numerical result for the percentage correction to the partial width T{H — )■ 77) has been presented 
in Ref. [38J in the minimal complex-mass (MCM) setup and is shown in Fig. [2H as well as the extension to 
the full complex-mass (CM) scheme. The QCD corrections (dotted curve) as well as the CM electroweak 
corrections (solid-line) are shown separately. If it is true that below the WW threshold (5qcd and ^ew almost 
compensate and lead to a small total correction (dash-dot line), the new results tell us that, above the WW 
threshold, both contributions are positive and lead to a sizable overall effect of approximately 4%. In the 
considered Higgs-boson mass range the total correction varies between — 1% < (5tot < 4%. 

To summarize, the electroweak correction around the WW threshold shown in Fig. [21] (dash-line) has 
been produced by minimally modifying the two-loop amplitude through the usage of a complex mass. 
We have extended our calculation in order to treat the whole amplitude completely with complex W and Z 
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Figure 21: The two-loop electroweak percentage corrections for the decay width T{H gg) and the total partonic 
cross section o{gg — > H). The solid line denotes the total electroweak correction, including also top quarks (first + 
second + third generation) . The vertical dash-lines indicate the location of electroweak thresholds. 
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Figure 22: The two-loop electroweak and QCD percentage corrections for the J/ — > 77 decay. 



masses. This allows us to study the remaining cusp and the influence of the complex mass on the electroweak 
correction, which is shown through the solid line. 

We have compared our numerical results in the region below the WW threshold with those of Ref. [28 . 
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Eqs.(8) and (11) of Ref. [28^ contain typosH ; once we correct, very good agreement between our calculations, 
up to 130 — 140 GeV, is found. To be more precise, we split the total contribution into light fermions and 
top (namely diagrams with at least one top line); for the light-fermion contribution the agreement is always 
very good because the result of Ref. [57] is exact, whereas, for the top part, the comparison gives a check 
of how good is the expansion of Ref. [28] . We confirm what is expected, for light Higgs-boson masses (say, 
around 115 GeV) the agreement between our results is again very good and starts deteriorating for heavier 
values (say, starting from 140 GeV). In Fig. [23] we have summarized our findings for i5ew in gluon-gluon 




Figure 23: The two-loop electroweak percentage corrections for H gg showing a comparison: total in CM scheme, 
light fermions only in RM scheme, light fermions only obtained from Ref. [27] (see main text for the error bars). 



fusion; here we also include the light-fermion part of the corrections compared with the ones of Ref. [27] . 
The result of Ref. [27j is known completely in terms of harmonic poly-logarithms. For simplicity we used 
the tool EasyNData [55] to read out the result from Fig. 2 of the first paper of Ref. , which is shown in 
Fig. [23] as a dash-line. The error bars originate from our estimation of how good we can read out the data 
from the plot. Once again, the relevant news are the behavior above WW threshold induced by top quarks 
and the around-threshold behavior in the CM scheme. 

Differences become even more striking when we consider the total corrections in a narrower range in the 
Higgs-boson mass, as done in Fig. 1241 Here we concentrate our attention around the region containing both 
WW and ZZ thresholds and compare our final result (CM scheme) with the one obtained from Tab. 2 of 
Ref. [28]. 

The very good agreement when we use the MCM scheme and the subsequent extension to CM scheme 
make us confident of the goodness of the result. 



9 Conclusions 

In this paper we have provided all technical details for a stand-alone numerical calculation of the full 
two-loop corrections to the decay widths ^ 77 and H gg, where H stands for the Standard Model 
Higgs boson. The techniques introduced in this context, however, have a much wider range of application, 

^Thanks to G. Degrassi and F. Maltoni for prompt confirmation. 
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Figure 24: An analogous plot to Fig. [23] for a narrower range of the Higgs-boson mass. Here we compare our results 
in the CM and MCM schemes with the results (below 155 GeV) of Ref. g8j 



i.e. they are general enough to be used for all kinematical configurations of 1 ^ 2 processes at the two-loop 
level. 

The generation of diagrams for any given process has always been performed with the FORM [39j pro- 
gram QraphShot [40] which performs simplifications and symbolic manipulations of the loop integrals, like 
performing traces, removing reducible scalar products and symmetrizing the integrals. 

The amplitudes contain both QCD and electroweak contributions and our strategy has been to present 
complete calculations where all enhanced terms are extracted analytically and their cancellation shown, 
whenever it occurs; the remainder of the amplitude is a multi-dimensional integral in parametric space with 
an integrand which, after suitable algebraic manipulations, is transformed into a smooth function therefore 
avoiding well-known problems of numerical instabilities, even in those regions of parametric space which are 
notoriously plagued by normal threshold singularities; the latter are always cured by the introduction of 
complex masses, solving at the same time a conceptual problem, internal unstable particles, and a practical 
one, avoiding numerical instabilities associated with the behavior around singularities lying on the real axis 
of the external Mandelstam invariant. 

The strategy for the rest of the calculations consists in regularizing the coUinear singularities and in 
extracting the singular pieces from all singular terms. After checking - analytically - that singular parts 
cancel in the total (if it applies), one can safely take rid of the regularization parameter and include all 
collinear-free remainders into the total amplitude. 

As far as numerical results are concerned we can safely state that the methods developed in this paper 
produce extremely accurate predictions for any value of the Higgs-boson mass, including the full dependence 
on the W-, Z- and Higgs-boson masses and on the top-quark mass. A consistent and gauge-invariant 
treatment of unstable particles allows to produce precise results around the physical thresholds: both a 
minimal and a complete version of this scheme are studied. It is found that the use of the complex mass 
scheme smoothens the threshold singularities and that a complete implementation of this scheme is needed 
in order to get the NLO electroweak corrections under control. 

When applied to the gluon fusion channel, pp^gg + X^H + X, our results show that the electroweak 
scaling factor which should multiply the QCD-corrected cross-section is between —4% and +6% in the range 
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100 GeV < < 500 GeV, without incongruent effects around the physical electroweak thresholds. Finally 
we observe that around the tt threshold the top-quark contribution to H ^ gg leads to a sizable effect of 
about -4%. 
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A Self-energies, vertices and tadpoles 

In this section we collect our conventions for the diagrams involved in the paper. 
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Figure 25: The one-loop self-energy and vertex. / is a generic polynomial in the loop momentum q. The dimension 
of the space-time is n = 4 — e and ^ is the renormalization scale. 
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Figure 26: The irreducible two- loop self-energies diagrams. / is a generic polynomial in the loop momenta qi and q2. 
The dimension of the space-time is n = 4 — e and fi is the renormalization scale. 
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Figure 27: The irreducible two-loop vertex diagrams. / is a generic polynomial in the loop momenta qi and 52- The 
dimension of the space-time is n = 4 — e and /i is the renormalization scale. 



B Properties of projectors 

In this appendix we briefly summarize a general approach based on the work of Ref . |52| . AmpHtudes 
for two-loop 1—^2 processes are decomposed into form factors which have to be extracted with proper 
projection operators. Let us consider tensor, one-loop, A^-point functions in n dimensions {N < 5, n = 4 — e) 

Si:,," = ^ f rqj^^—^^ {i) = {q+Vl+ ■■■ +V^f + ml (234) 
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Figure 28: The irreducible one- and two-loop vacuum diagrams. The dimension of the space-time is n = 4 — e and 
is the renormalization scale. 



where po — 0. We select — 1 independent vectors, n, i = 1, . . . ^ N — 1 and introduce the following 
notations: 



G,;; = G„ ; = 2 r, • rj, 5" = det 



For definiteness we will choose — pi. Next we introduce 

1 



n- 7V + 1 
They satisfy the following properties: 



^^iU 2 T ^ G ^ Ti/ 



Let us consider now the action of these projectors on tensor integrals: we consider first 
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with xq — 1 and x„ = 0. Introducing x ({2;}) = Af ^ — P ■ P, we obtain 

27 2 

When the integral is projected with D we have 
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For scalar diagrams we define the fundamental parametric representation, 

X = X -I- 2 x + L„_i = (x - X)* (x - X) 

Furthermore, we define new scalar objects, 
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/ 2\2-n/2 

where (2) = B, (3) = C etc. We obtain 

^i^i.R'jv ^u'n — S'„jv(0) + - Sn+2 N{h j)- 

Next we consider integrals with three momenta in the numerator 
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where fully symmetrized tensors have been introduced. There are two kind of projections, 
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With four momenta in the numerator we have 
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where fully symmetrized tensors have been introduced. Projections to be considered are: 
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K^I^RL Si:7 = - SnAm - G;^, 5„+2„(fc) - G-]^, 5„+2«(j) - G-]^, 5„+2„(z). (248) 
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C Non-abelian soft and collinear diagrams in H gg 

In this appendix we present a short summary of our procedure for extracting soft/coUinear divergent 
parts of the two-loop diagrams that contribute to the process H —f gg. Mass-singular, virtual, configurations 
appear when a) a massless, internal, line is connected to two on-shell, external lines (soft); b) a massless, 
external, line is connected to two massless, internal lines {collinear). Furthermore, we have a soft singularity 
whenever we attach an extra massless line to an on-shell external line; if the latter is also massless an 
additional collinear divergence arises. These QCD-like configurations have been extensively discussed in the 
literature, e.g. see Ref. [5], and in this section we briefly illustrate our approach to the problem. 

It is worth noting that in the QCD sector of the corrections to H ^ gg we select the regulator of 
collinear divergences to be the space-time dimension and not the masses of the light quarks. On the contrary 
we select the latter in all situations where collinear singularities cancel in the total. In the following we 
adopt the convention that a letter i inside a one- loop diagram denotes the scalar product q ■ Pi, where q is 
the corresponding loop momentum. 

In the soft /collinear decomposition of a two- loop diagram we could highlight the nature of the singularity 
by introducing factors F^"^ = and write a soft/coUinear decomposition for an arbitrary two-loop vertex, 

^ib= E K^■,^P-1iMl), (252) 

1=0,2 
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where the superscript / denotes the family whereas the subscript lab denotes the rank of the tensor (0 for 
scalar, etc.) as given in Ref. 49J. For those integrals where also ultraviolet divergences are present we need 
to extend the decomposition of Eq. (I252p to include higher orders of e poles by introducing 



Fl.i.) .4 + ^1^ + \aU^). (253) 



Alternatively, we can factorize UV and IR/coUinear divergences. Suppose that we are considering a two-loop 
diagram with a single UV pole; remembering that each loop can contribute at most one soft and coUinear 

term the highest possible pole in our integral is Consider the case where the highest pole is 

then we first give the UV decomposition and further decompose the coefficients with new IR factors: 

1=0,3 4=0,1 

Ki^I = V^^,,f'T{mII V^:;^ = E V^^.F^riMl), (255) 

1=0,2 



F^'-'"'{x)=Frix), FJ, (x) = -- + -A^„,(x). (256) 
In the following we give a sample of our results; the full list can be found in |92j. 



m 



H 



V 





H 



Fieure 29: The V'", V^'^ and V" non-abelian soft/collinear configurations. Solid lines stand for a common mass 
m. 

For the y^-" configuration we introduce a quadratic form, x{x) — x{x — 1) + m? /M"^, where m is the 
mass in the triangle of Fig. [211 and obtain 



y"2'" = 4 / dx-, 
Jo X 

y-{ = 2 / - j fdy In (x-y(:r- jy)) - ^ In - ^^^^ | + x 



X 



\n{x-x{l-x))-'i\nx+iTT k(257) 



For the V"^'^ configuration, we introduce a new quadratic form, ^{x,y) = x{^) + 2/^™^/-^^7 and derive 
the following result: 



l^_^2'' = 0, r^i''= I dx I dy 



ln(l - y) 



(258) 



Integrating over y we obtain 
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(259) 
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with y± — 1± I (ra ^ — x ) • 

Finally, for the non-planar V" configuration of Fig.[2Sl we obtain the following result: 



Jq X \ / 



Ml 



(260) 



and carry on the x integration, 
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1,2 
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(261) 



with x± = 1/2 ± 1/2-^1 — 4m2/il/2 . Here Li2 (z) and S'i.2(z) are Nielsen poly- logarithms. 



D Loop integrals with complex masses and momenta 

Consider the logarithm of a complex number 2 = z„ + i Zj. Let z = 2^ " * 0, we define 

ln(2) (2; f) =lnz-2z^0(-2«), (262) 

which satisfies 

lim ln(2) (2 ; 2) = In (2) . (263) 

When computing the amplitude in the CM scheme we will encounter, after reduction, one-loop two-point 
functions where both masses and the external invariant are complex. Let 

X{x) ^spx^ + {Ml - Ml - sp) x + Ml, (264) 

= Af 2 - i r M, Ml ^ml -i -fi rui . (265) 
The correct definition of the _Bo sfunction is as follows: 

1 

So = - - / dx ln(2) (x ; X) , X = X - * 0, (266) 

e Jo r,7i=o 

with the only restriction that x(a;) does not cross the negative real axis for x e [0, 1] (in which case the 
definition itself of a two-point function becomes dubious). The best way of understanding the second Riemann 
sheet is to consider complex and zero masses; here x = x{\ — x), therefore Re x(2;) < but Im x{'^) > 
which is the opposite of the — i prescription. Thus In gives the wrong answer and ln(2) the correct one 
based on the fact that we require a continuous limit widths — *■ 0. Therefore, to compute Bq one should start 
with the definition of the corresponding function with masses and real and continue to complex masses 
and p^ on the correct Riemann sheet, i.e. with ln(2) instead of In. Note that in this function y ± i is always 
treated as y ±i. Of course, one could use directly the logarithm of a complex number to compute In a;. 
When only the internal masses are complex there is no problem at all; indeed we find 

Im X = —71 mi (1 — x) — 72 m2 a; < 0, Im x = —0 < 0. (267) 

This fact is also true for arbitrary two-loop diagrams as long as only internal masses are continued into the 
complex plane. There are special cases of complex p^ , complex masses. In the actual calculation we have at 
most two scales: consider this situation where we are free to set p^ = —Ml and obtain 

x{x)=Mlx^ + Ml{l-x), (268) 

which means that Re x{^) > ^ ^ [O7 1] and, therefore, we can use standard results, i.e. ln(2) In. 
Furthermore, when Mi = we obtain 

x{.x)=p\l-x) + Ml (269) 

If p2 = —M'^ + i F M then Re x{^) < requires M > 1112 which is never satisfied in our case (M^ < Mt, M^). 
Also here ln(2) — > In. Therefore, as far as the calculation reported in this paper is concerned, the only relevant 
case is complex p^ and zero internal masses. 
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